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abstract 

The  analysis  presents  a practical  method  for  evaluating  the  added 
mass  and  damping  coefficients  of  a heaving  surface  effect  ship  in  uniform 
translation. 

The  theoretical  added  mass  and  damping  coefficients  and  the  heave 
response  show  fair  agreement  with  the  corresponding  experimental  values. 

Comparisons  of  the  coupled  aero'hydrodynamic  and  uncoupled  analytical 
results  with  the  experimental  data  prove  that  the  uncoupled  theory,  dominant 
for  a long  time,  that  neglects  the  free  surface  effects  is  an  oversimplified 
procedure. 

The  analysis  also  provides  means  of  estimating  the  wave  elevation  of  the 
free  surface,  escape  area  at  the  stern  and  the  volume  which  are  induced 
by  a heaving  surface-effect  ship  in  uniform  translation  in  otherwise  calm 
water.  Computational  procedures  have  been  programmed  in  FORTRAN  iV  language 
and  adapted  to  the  POP-IO  high-speed  digital  computer. 


KEYWORDS 

Surface  Effect  Ship 
Water  Waves 

Aero-hydrodynamic  Forces 
Bubble  pressure 
Added  Mass 
Damping  Factor 
Heave  Response 
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NOMENCLATURE 

A escape  area  or  area  of  nozzle  throat 

A^  escape  area  at  the  equilibrium 

A cushion  area 

c 

A^  variation  of  the  seal  opening  area 

a half  length  of  the  pressure  patch  or  wave  amplitude 

b half  width  of  the  pressure  patch 

C_  non-dimensional  added  mass  coefficient 

O 

Cj  non-dimensional  damping  coefficient 

Ff,  Froude  number 

g gravitational  acceleration 

h cushion  height 

i ^ , or  the  suffix  indicating  the  imaginary  component 

j integer 

K nozzle  flow  coefficient  or  constant 

kj  wave  number 

L nominal  ship  length 

I integer 

H"  added  mass 

m integer 

N damping  factor 

n integer 

P pressure 

p atmospheric  pressure 

P|^  bubble  pressure 

p^  ambient  pressure 
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c; 


K 


Q discharge 

discharge  at  equMibrIun 

r distance  or  the  suffix  Indicating  the  real  component 

s distance 

m 

t time 

U ship  speed 

V the  wave- induced  volume 

V non-dimensional  coefficient  of  the  wave- induced  volume 

w suffix  indicating  the  water 

x,y,z  coordinates 
z heave  displacement 

a non-dimensional  coefficient  of  the  escape  area 

V ratio  of  specific  heats 

c a parameter 

9 the  velocity  potential 

\ scale  ratio 

n pressure 

a circular  frequency 

T frequency,  oU/g 

p,p^  density  of  water 

p density  of  air 

C wave  elevation 
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INTRODUCTION 

The  analysis  is  directed  at  the  problem  of  finding  the  effect  of  the 
presence  of  the  water  wave  upon  the  bubble  (or  cushion)  pressure  in  the 
plenum  of  a surface-effect  ship  as  it  is  forced  In  simple  harmonic  heaving 
motion  while  translating  at  constant  speed  over  otherwise  calm  water. 

The  surface-effect  ship  undergoing  the  above-mentioned  motion  is  a 
mult {-parameter  system  influenced  by  the  coupled  aero-  and  hydrodynamic 
actions  of  water  and  air  in  the  plenum  which  is  under  a continuous  supply 
of  air  through  a system  of  fans  and  leakage  through  the  peripheral  gap. 

The  analysis  when  linearized  yields  the  aero-hydrodynamic  forces  due 
tc  the  forced  heaving  motion  and  thus  provides  the  corresponding  added  mass 
and  damping  force  coefficients  in  waves  which  is  vital  information  for  the 
evaluation  of  particular  design  characteristics  of  a ship  at  the  desired 
speed  and  sea  state. 

The  theory  has  been  developed  by  referring  to  the  works  by  Kaplan  and 
12  5 

Davis,  Doctors,  and  Breslin  and  Hires.  The  relevant  experimental  litera- 
ture referenced  here  are  the  investigations  due  to  Van  Den  Brug  and  Van 

Staveren,^  Hagnuson  and  Wolff, ^ Fein  and  Murray,^  Moran,  Fein  and  Ricci, ^ 

8 9 

and  Fridsma.  ’ 

Kaplan  and  Devis^  assume  that  the  free  surface  behaves  as  a rigid 

boundary  at  constant  pressure  whereas  Doctors  and  Breslin,  et  al,  assume 

that  the  deformation  of  water  surface  participates  in  the  generation  of  the 

bubble  pressure  in  conjunction  with  the  actions  of  seals  and  fans,  etc.  The 

2 T 

second  cardinal  assumption  in  the  latter  works  ’ is  that  the  deformation  of 
the  wave  surface  under  the  oscillatory  rectangular  pressure  patch  in  uniform 
translation  be  used  to  display  the  way  in  which  the  motion  of  the  water  surface 
participates  in  the  determination  of  the  pressure  variations  in  the  plenum  air. 

Doctors  deals  with  the  three-dimensional  rectangular  pressure  patch 
whereas  Breslin  and  Hires  idealize  the  problem  to  two-dimensional  flow  condi- 
tions. 

^Superior  numbers  in  text  matter  refer  to  similarly  numbered  references  listed 
at  the  end  of  this  report. 
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The  analysis  has  been  developod  Independently  from  that  of  Doctors, 
a|thou9h  along  the  same  lines,  leading  to  expressions  for  the  evaluation  of 
the  wave  elevof  i,  t sc.ipe  area  and  volume,  quite  different  from  those  re- 
ported in  Refeit-Mv^t  2.  In  addition,  a nun^erical  approach  has  been  developed 
valid  for  the  entire  range  of  frequency  of  practical  interest,  from  very  low 
to  considerably  high,  in  contrast  to  Doctors*  numerical  results  which  are 
restricted  to  the  very  low  frequency  regime. 

In  the  dynamic  analysis  of  the  airflow  phenomena,  pressure  and  volume 
changes  are  assunted  as  by  others'”^  to  occur  sufficiently  rapidly  so  that  the 
adiabatic  law  qoverns  the  basic  thermodynamic  variations  in  the  air  cushion. 
The  cushion  pressure  is  assumed  to  be  spatially  uniform  which  is  approximately 
valid  according  to  Moran'^  and  Tsakonas  et  al.*'  Since  this  work  is  intended 
to  be  applied  to  the  high  length-beam-rat lO  sidewalled  surface  effect  ships, 
it  assumes  also  that  the  leakage  of  ai r takes  place  only  at  stern  seal.  The  seal 
response  dynamics  are  beyond  the  scope  of  the  present  investigation.  It 
assumes  a semi-rigid  planing  seal  which  has  been  the  subject  of  the  experi- 
mental study  of  the  XR-5  SES  Model. ^ 

Following  the  calculation  of  the  waves  due  to  an  oscillatory  rectangular 
pressure  patch  in  uniform  translation,  the  bubble  pressure  is  determined.  The 
linearized  equation  of  bubble  pressure  in  plenum  chamber  and  that  of  the  heav- 
ing displacement  are  both  derived  by  applying  the  perturbation  method  through 

’ the  time  derivative  of  the  adiabatic  law.  The  asymptotic  behavior  of  the 

i 

bubble  pressures  as  the  frequency  -*  0 and  » is  analyzed  and  discussed  with 
and  without  water  wave  effects. 

The  mot  ion- induced  aero-hydrodynamic  forces  are  calculated  from  the 
bubble  pressure  distribution  on  the  cushion  area  as  a function  of  Froude  num- 
ber and  of  reduced  frequency.  From  its  real  and  imaginary  components,  the 
added  mass  and  damping  coefficients  are  determined  respectively.  The 
asymptotic  behavior  of  the  added  mass  and  damping  as  the  frequency  -*  0 and  «• 

! are  evaluated  and  discussed  both  with  and  without  water  wave  effects. 

I The  scale  effect,  essential  Information  for  the  extrapolation  of  the 

I model  measurements  to  the  full-scale  condition,  is  also  the  subject  of  the 

I 12 

present  analysis.  The  scale  effect  determined  in  the  previous  study  in- 
j dicates  a great  dissimilitude  in  general.  The  scale  effect  ratio  of  the 

I 
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pressure  of  the  ship  to  the  model  as  calculated  here  shows  that  the  dis- 
similitude is  small  at  low  frequency  but  quite  large  at  higher  frequencies. 

The  present  analysis  does  not  include  the  sid>.'w^all  hydrodynamic  ef- 
fects since  experiments  on  the  XR-5  SES^’^  have  shown  slight  effect  of  the 
sidewalls  on  added  mass  and  damping. 

A series  of  numerical  computations  has  been  done  to  determine  the  water  ' 

waves,  bubble  pressure,  heave  added  mass  and  damping  coefficients,  scaie 
effect  on  bubble  pressure,  and  the  heave  responses  of  the  XR-5  SES  model  run- 
ning in  head  waves. 

The  waves,  escape  area  and  voiume  which  are  induced  by  the  oscillating 

pressure  at  uniform  speed  in  otherwise  calm  water  have  been  calculated  for 

2 

the  model  used  by  Doctors.  A computational  procedure  has  been  developed 
and  established  for  any  frequency  and  speed  without  limitation.  The  bubble 
pressure,  added  mass  and  damping  coefficients  have  been  computed  for  the 
XR-5  SES  model,  with  and  without  incorporating  the  water  elevation.  A 
remarkable  difference  has  been  found  between  the  two  groups. 

Experimental  values^  of  the  heave  added  mass  and  damping  coefficients 
have  been  compared  with  the  theoretical  results.  To  study  the  validity  of 
the  theoretical  added  mass  and  damping  coefficients,  the  heave  response  of 
the  model  has  been  calculated  on  the  assumption  that  the  pitch  is  restrained, 
using  the  theoretical  added  mass  and  damping  factors  with  the  experimentally^ 

I determined  static  stiffness  and  heave-exciting  forces. 

The  analytical  developments  at  limiting  cases  (i.e.,  high  and  low  fre- 
quencies) have  been  derived  and  discussed  in  detail  as  well  as  the  results 
of  the  comparison  of  predictions  with  measurements.  The  comparisons  of  the 
results  of  both  the  coupled  aero-hydrodynamic  analysis  and  the  uncoupled 
analysis  with  the  results  of  experiments  have  shown  strongly  that  omission 
of  the  water  waves  in  calculating  the  bubble  pressure  and  the  hydrodynamic 
forces  is  not  a valid  assumption. 

This  research  was  sponsored  by  the  David  Taylor  Naval  Ship  Research  and 
Development  Center,  General  Hydromechanics  Research  Program,  under  Contract 
N000I4-77-C-0061 , Davidson  Laboratory  Project  4480/015. 
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THE  velocity  potential 

The  hydrodynamic  boundary  value  problem  for  a pressure  patch  oscillating 
with  constant  amplitude  and  moving  on  the  surface  of  a deep,  incompressible, 
ideal  fluid  can  be  determined  through  the  velocity  potential  method. 

The  right-handed  Cartesian  coordinate  system  is  set  up  on  the  calm 
water  surface  with  its  vertical  z-axis  pointing  upwards.  See  Figure  la. 

13  14  , 

This  problem  is  formulated  as  an  unsteady  boundary  value  problem  ' in 
the  framework  of  the  linearized  conditions. 

The  pulsating  pressure  distribution  of  a rectangular  shape  on  the 
calm  water  surface, which  is  a simple  harmonic  in  time,  is  expressed  as 


p(x,y,0;t) 


n(x,y)e 


txl  ir  a,  lyl  ^ b 
Ixl  > a,  lyl  > b 


. t > 0 


The  deforma tion  of  water  due  to  the  above  pressure  distribution  is 
assumed  to  be  small  and  is  represented  by  the  velocity  potential  v(x,y,z;t) 
which  is  harmonic  in  space  below  the  free  surface;  satisfying  the  d.e., 

9+cp+cp-O,  z<0,t>0  (2) 

XX  yy  zz 

On  the  free  surface  the  velocity  potential  and  the  surface  elevation 
C ■ C(x,y;t)  satisfy  the  linearized  kinematical  and  dynamical  boundary  con- 
di tions, 

-Ct  + UCx  + «Pz  - 0 . ^-0.  (3) 


•£  + qC  + to  - 
p t 


z-0,  t>0 


Eliminating  Q from  Eqs.(3)  and  (4),  we  have  the  linearized  free  sur- 
face condition 

'^tt  ^‘^xx  ■ ^^xt"^  99j-^(Up^-Pj).  z-0  . t > 0 (5) 

In  addition  to  the  above,  «p  and  satisfy  the  condition  of  being 
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initially  at  rest  on  the  surface 

cp(x,y,0;0)  - «P,.(x,y,0;0)  - 0 (6) 

where  9^*0  is  valid  under  the  assumption  p(x,y,0;0)*0  [Eq.(1)],  and  on 
the  bottom 

- 0 , z - t > 0 (7) 

In  the  far  field  we  suppose  that  cp  and  its  first  and  second  deriva- 
tives tend  to  zero  for  any  given  time;  in  fact,  they  tend  to  zero  in  such 
a way  that  Fourier  transforms  exist. 

Now  we  seek  the  form  of  cp  which  satisfies  the  above  conditions  by 
employing  Fourier  transform  technique.  Define  the  double  Fourier  trans- 
forms 

cp(w,u,z;t)  = J cp(x,y,z;t)e"'^'^^‘^“^^dxdy  (8) 

**00 

and  inversion 

cp(x,y,z;t)  = J 9(w,u,z;t)e' ^'^^■^“^^dwdu  (9) 

•eo  -00 

Here  w and  u are  the  longitudinal  and  transverse  wave  numbers  and  are 
related  to  the  circular  wave  number  k and  the  wave  angle  0 by 

w “ k co-»0 

u “ k s i n0  (10) 

The  Fourier  transform  (8)  applied  to  Eq.(2)  yields 

“ 0 (11) 

The  general  solution  of  the  above  differential  equation  is 

9 * A(w,u;t)e’*'*^+  B{w,u;t)e"*’^*^  (12) 

where  A and  B are  arbitrary  constants.  By  making  use  of  the  Fourier 
transform  of  the  bottom  condition  (7),  we  obtain  B=0  . Hence  the  above 
is  reduced  to 

cp  ■ A(w,u;t)e**^*^  ( 13) 
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To  determine  A,  take  the  Tourier  transform  of  the  linearized  free  surface 
condition  (5), 

<ftt  “ - iZUwcp^  + g(p^  - i(Uw-o)  ~ e*°*  (14) 

where 

n (w.u)  - ~ J J'  n(x,y)e"’ 

dT)J  n (5.T))e"'‘^‘"5+uTl)j^ 

-b  -a 

Substituting  $ from  Eq.(13)  in  (14)  yields 

A^^-U  w*A  - 12UwA^  +g|k|A-  i(Uw-cj)  ~e  (16) 

which  by  means  of  the  Laplace  transform  and  utilizing  the  initial  condition 


A(w,u;0)  " A^(w,u;0)  - 0 

yields  the  following  solution  in  the  transform  plane 


L(A). 


( Uw-a) n 


Zpvgikl  '’"‘’l  ‘’“'*2 


wl  th 

qj  - i (Uw  + Vglk'l  ) 
q^  - i (Uw  - Vg|k|  ) 


Inversion  of  the  above  is  given  by  the  convolution  integral, 


A-:tl  2p  vg|k|  o 

Equation  (17)  after  integration  is  substituted  in  (13)  and  we  have 

ip  - jklz[l-e'^^'*‘^^^'~^^^J 

^±1  2p  ^/g]^^^  Uw  + Jt  -Tgirr  " o 


(17) 


(18) 
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The  inverse  transform  of  (18)  is  given  by 


^ I J'  (5,^1)  J J 

P •h  wA  •<»  •a>  . /i^  iL 


? 1 |h|i+i  (wx_+uy„) 

P ^-n  /•  'n\  P P ' ' ' O 'o' 


P -b 


-v/gl^l 


(uw-o) 


, _i(Uw  + Vgikl  -0)t  , i(Uw -v/gTkr“a)t 

I “•  I "0 


Uw  + Vg  I k I “ o 


uw  - 7^ 


~) 


dwdu  (19) 


where  x ■ x - C 
o ’ 

Vo  - y - Ti 


Now  the  variables  (w,u)  are  transformed  into  (k,0)  by  the  relation  (10) 
Replacing  the  elementary  area  dwdu  with 


dwdu  ■ kdOdk 


for  the  intervals  0 s k s « , -n  s.  0 s n,  we  have 


cp  - e 


iot  i 


J dll  J'  d^n  (?,T1)  J d0  J dkJE  (ukcosO-o) 
-b  -a  -TT  o * ^ 


(10') 


kz+i k(x  cos0+y  s i n0) 


1 _ pi  (Ukcos0+ Vgk  -a)  t | _ pi(Ukcos0  -o)t 


Ukcos0  + VglT  -a 


Ukcos0  - “ o 


To  fold  up  the  integration  interval  (-n,n)  to  (0,  y)  , the  integral  J 
is  subdivided  into  four  quadrants, 

n -n/2  o tt/2  n 

P - P ^ P,  * P * P, 

_n  -n  -n/2  o n/2 

n/2  -n/2  n/2 

and  each  Integral  is  transformed  into  J . For  instance,  J -•  J*  is 

o -no 

achieved  by  setting  0«Y-n  , where  Y is  a new  variable.  It  follows  that 
cot0  ■ -cosY,  sln0  ■ -sInY  . 

o n 

By  replacing  Y with  6 and  applying  a similar  procedure  to  J and  J , 
we  have  the  following  folded  form, 
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cp  ■ e 


lot  I 


tt/2 


J dTl  J dsn  (5.T1)  J de  J dk  J| 


8n*  p »b  -a 


o o 

iu>i  t 


|(UkcosO-o) 


1-e 


iu^t 


(Ug 


« _x  / -ikr*  -IkraW  l-e 
+(Ukcose+o)  (e  ^+e  )\-  — — + 


ioiat 


UJ4 


■)} 


where 

■"i  “ (x-g)cosW  + (y-'n)sin0 
ra  ■ (x-g)cos0  - (y-Tl)sin0 

and 

'‘*1  'I 

I ■ Ukcos0  ± - o 

u)a 

“ Ukcos0  ± + a 

(i)4 


(22) 


(23) 


(24) 
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STEADY  STATE  SOLUTION 

In  Eq.(22),  cp  is  the  solution  of  unsteady  conditions  which  have 
been  stated  at  the  outset.  The  next  step  is  to  derive  the  steady  state 
solution  from  Eq.(22)  by  eliminating  the  time-dependent  terms.  This  is 
achieved  in  the  following  manner;  1)  determine  the  wave  numbers  k from 

I 

'".“0  , and  2)  take  the  proper  indenting  paths  in  the  neighborhood  of  k. 
in  the  k-plane  in  order  to  have  the  exponentials  vanish  in  time  as  t -•  «». 
The  integration  paths  determined  as  such  are  therefore  not  arbitrary  but 
strictly  bound  to  derivation  of  the  steady  solution  and  to  the  roots  k,. 

Solving  for  the  zeros  of  in  Eq.(24),  namely  k,,  we  have 

v2 


or 


^1 

k„ 


(ii 


1 + 


4Uocos9 


) 


k[f  cos  9 


(25a) 


(i  1 - hiESPii  'f' 


4u^cos  9 


(26a) 


o J 2 


k^sec^9  [1+2t  cos9  + V1+4tcos9  ] 


(25b) 


J j k^sec^9  [1-2t  cos  9 + V1-4tcos9  ] 


(26b) 


with 


i- 

If’ 


9 


(26c) 


It  is  seen  from  (25a)  that  kj  and  k^  are  positive  real  and  k^*  k2 
for  interval  0 ^ 9 g ^ , it  is  seen  from  (26a)  that  k^  and  k^^ 
are  positive  real  and  if 

j . 4Uqcos9  > Q 
9 
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and  k-  and  are  complex  with  |mk-<0  and  imk.>0  (i.e.,  k.“  complex  conjugate 
of  k J - 


If 


1 - < 0 

9 

In  other  words  k^  and  k^  are  positive  real  in 

0 < e < 5 

c « =2 

(27) 

and  complex  in  the  interval 

0 < 0 < 0 
“ c 

(28) 

where 

®c  ■ <i;u5> 

(29) 

Now  we  consider  the  integrals  with  respect  to  k in  the  k-plane. 

Since  Taylor  expansions  of  k and  Vgk  in  the  neighborhood  of  k., 
up  to  the  first  order,  are 

k - k.  + (k-k.)  + 0(k®) 


<ygk  z Vgk.  + (k-k.) 

2v/ir 


(i=J.2,3,4) 


(30) 


we  obtain  the  following  expressions  in  the  neighborhood  of  k. : 

i (Uk,cos9^^/gi^^-a)  + i (k-k.)  (Ucos9  + — ) 

2 ^ 

i (Uk2Cos&<-'/gi^-a)  + i (k-k2)  (UcosG  + — ) 


i (Ukcos0+v§l<-a)  - | 


2 ^/i^ 


I (UkcosG- v^-o)  - 


i (Uk^cos9-^/gk]  -<j)  + l (k-ki  ) (UcosG  - — ) 

2 


1. 


i (Uk2COsG-V^-a)  + i (k-k2)  (UcosG  - — 2^)  (31) 


2 ^/k^ 
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i(Ukcos0-»-  Vgk  +o)  - 


i(Uk,co50+  ^/giTT  +o)  + i (k-k.)  (Ucos8+  — ) 

^ ^ ^ 2-yi^ 


I (Uk^^co50+  ^/gk]^  +o)+i  (k-k^^)  (Ucos6+ 


2 


l(Ukcos0-  +o)  - 


i (Uk-cos0-  <ygk,  +o)  + i (k-k.)  (Ucos0 ) 

^ ^ ^ 2VkJ 

i (Uk.  CO50-  -y/gk.  +o)  + i (k-k.  ) (Ucos0 ) (32) 

^ 2Vi^ 


From  Eqs.(2Sa)  and  (26a),  we  have  the  positive  real  root  k.  for  the 
interval  0 ^ 0 ^ j 


7g  ^-1  +\|  I + 


♦UOCOS0 


ZUcosO 


Vg  ^1  + >]  1 + 


♦UOCOS0 


2UCOS0 


and  for  the  interval  0 S ® = T 

C ” fc 


^/gVl  -Nl- 


♦UOCOS0 


ZUcosO 


*^1  +“^1  - 


+UOCOS0' 


2UCOS0 


Employing  Eqs.(33)  and  (34),  we  have 
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l(k-kp 


./g  N 1 + 


l(Ukcos0+  -o) 


+ i2  M^gk^ 


4UOCOS0 


2 ^ 


i(Ukcos0-  Vgk  -o) 


•i2  ^/ik[ 


'(k-k.) 


^ >|  V .H 
2 ^ 


I (Ukcos0+  Vgk  +a)  "* 


i2  >/gi^ 


i2  Vikj; 


i (Ukcos0-  /v/gic  +o)~ 


1 

-l(k-k ) ^ 

^ 2 ^/i^ 


+i(k-k^) 


4UOCOS0 


^/g  N 1 - 


4uacos0 


zVki; 


Equation  (35)  shows  that  for  the  interval  0 ^ 0 ^ Y the  path  of  integra- 
tion of  the  terms  containing  e'**^^  and  e'***®^  in  Eq.(22)  in  the  complex 
k-plane  should  be  deformed  above  the  real  axis  in  the  neighborhood  of  both 
k|  and  k^  in  order  to  have  the  exponentials  in  t vanish  as  t -•  « . Let 
this  deformed  path  be  denoted  by  L p (See  Fig.  1b.)  Equation  (36)  shows  that 
for  the  interval  0^.  i 0 < ^ • th®  path  of  integration  containing  the  factors 

e and  e in  Eq.(22)  in  the  complex  k-plane  should  be  deformed 

above  the  real  axis  in  the  neighborhood  of  k^  and  below  the  real  axis  in 
the  neighborhood  of  kj^  In  order  to  have  the  exponentials  in  t vanish  as 
t -•  «.  Let  this  path  be  denoted  by  L2  • Equation  (24)  shows  that  for  the 
complex  wave  numbers  k^  and  kj^  in  the  interval  0 it  0 < 0^,  the  integration 
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path  requires  no  deformation  In  order  to  have  the  exponentials  in  t 
vanish  as  t - « ; for  the  path  along  the  real  axis  the  exponentials  vanish 
as  t - Let  this  be  denoted  by  Lj  . (See  Figure  lb.) 

From  Eqs,  (2A),(25a)  and  (26a),  we  have  the  following  relations, 


- i TiT;  u=co,»e  r-U  ^ 

2 'u)3  \k-kj  k-kj^ 


k-k,  ) 


(37) 


Then  substituting  the  above  relations  in  Eq.(22)  and  omitting  the  time 


factor  e'^*  we  have, 


b a 


Tl/2 


{-  J d0sec“e  J dk  F,(e,k) 


o 

L, 


TT/2 


+ r dGsec^e  J’  dk  F2(9.k)  + J dOsec^e  J dk  F,(e,k)} 

Q f\  r\  ^ ^ 


(38) 


where 


F,(e,k)  - (ukco59-o)(e''^*'*‘''‘iUe‘'^*'*‘*'‘»^)  — - Jli.  J | 

^l”^2'^”^l  k-k^'  I 


The  above  cp  is  the  ultimate  form  of  the  velocity  potential  for  the 
steady  state  oscillation. 
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SURFACE  ELEVATION 


The  surface  elevation  at  a point  on  the  calm  water  surface,  z»0. 

Is  evaluated  by  assuming  uniform  pressure  distribution  II  - 11(5,11)  of 
a rectangular  shape,  omitting  the  time  factor  e**^^  , 

^ ~ 'P  . *-0  (40) 

where  the  first  term  is  the  static  elevation  while  the  remainder  represents 
the  elevation  due  to  the  dynamic  effect. 

Let  the  velocity  potential  cp  in  Eq.(38)  be  written  in  the  follow- 
ing form,  omitting  the  factor  e'^^  — S— — , 

4n*pu^ 

Ti/2  . k,  k 

»--J  d6,ec*e  |(8)  -^9, 


n/2  k,  k. 

* i 


»C  k k 

d9s.c=9  1^jpl-93  , (8)  -5^93,2(8); 


wi  th 


b a 


<Pl  j(8)*  r dll  J d5  J dk(Ukcos0-o) 

’ -D  -a  o 

Ll 

b a “ 

2^®)"  I J *^5  J dk(Ukcos9-o) 

’ -o  -a  o 


k(2+irj)  k(z+ir2) 


+e 


k-k. 


k(z+irj)  k(z+ir2) 


+e 


k-k„ 


k(z-ir|)  k(z-ir2) 


b a • 

*^2  1^®^*  ^ J J dk(UkcosO+o) 

’ -b  -a  o k-k 


+e 


(41) 


b a » k(z-ir,)  k(z-ir2) 


«p2  2^®)"  / J *^5  I dk(Ukcos0+o) 


+c 


-b  -a  o 


k-k, 


(42) 

CCont 'd] 
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b a « k(z-lr,)  (.(z-ir^) 

V,  ,(9)-  J dTl  J d§  J dk(Ukcos04o)  S 

* -b  -a  o 


u a 

tpj  2^®^”  J I J'  ‘**<(Ukcose+o) 

’ “b  -a  o 


k-k. 


k(z-ir,)  k(z-ir,) 


+e 


k«k, 


('♦2) 


where 

fj  « (x-5)cos©  + (y-T])s(n© 
r2  ■ (x-5)cos6  - (y-TDsinG 

The  first  subscript  In  cpjj  denotes  the  path  of  Integration  as  Identified 

In  Figure  2 and  the  second  denotes  the  singular  point  along  the  given  line 

from  left  to  right.  To  evaluate  the  Integrals  systematically,  rearrange  the 

above  In  the  following  form:  a 

b a « •‘Vz+Tj; 

1 , ,(x.y.z;e,k,)  - J'  dT]  J d^  J dk(ukcose-a)  P 


-b  -a  o - 
-yv. 


b a » 

l|  |(x,y,2;0,kj)  » J dll  J d5  J dk(Ukcos0-a) 


-b  -a 


k{z+lr2) 

k-k. 


J,  ,(x.y.2;0,k,)  I,  , 

J,  ,(x,y,z;0,k,)  - ^ I 


(‘43a) 


b a » Mz+ir,) 

' 1 “ J dll  J d5  J dk(Ukcos0-a)  - — 

-b  -a  o 2 

_ru 


k(z+lr2) 


j I ® 1 1<2^ 

- r dl)  J d5  J dk(Ukcos0-a)  r— 

-b 

-a  o '2 

-ru 

9 

'l,2 

] 2(x.y.z;0,k2) 

d 

" 37 

'l.2 
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b a ™ 

I-  ,(x,y,z;e,k-)  “ J*  J J dk(Ukcos0+a) 

^ _b  -a  o 

_ru 


b a ® 

l'  ,(x,y,z;0,k-)  - J*  <m  J d?  J dk(Ukcose-Kr) 

^ _b  -a  o 

_ru 


'^2  I y ^ ® * ^3^  " 53r  ^2  1 

>2,1 


k(z-irj) 

a 

k-k. 


^(z-ir^) 

k-k. 


k(z-i Tj) 


b a « 

1^  2(j<.y.z;0.'^4)  = I I J dk(Ukcose+a) 

•b  ^ ^ 


b a « • 


k(z-ir2) 


'2  ^ I '‘S  J dk(Ukcos0+a) 


"vr 


J2  “ dx  '2,2 

jJ  2(x.y.z;6»kj^)  * 1;^  '2,2 


b a OD 


k(z-irj) 


I3  j(x,y,z;0,k^)  * J dll  J d?  J dk(ukcos0+a)  — 


-b  "3  ° 


^(z-ir^) 


b a “ 

I3  ,(x,y,z;0,kj)  » J dll  J d^  J dk(Ukcos0+a)  - 


J3  ,(x.y.z;0.k3)  I3  , 


J3  ,(x,y,z;0.k3)  = I3  , 
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b a k(2-ir,) 

. f dTl  J d5  ; dk(Ukcos0-Kj)  ^-j-r 

-b  “fl  o 4 


b a 'o 


■ J dTj  J dj  j dk(Ukcos6-fa) 


-b  -a  o 

a , 

*3.2 


Jj  2^**y'^*®*'^4^  - 1;^  I3  2 


(43f) 


The  Integrations  in  Equations  (43a)  through  (43f)  are  given  in 
Appendices  A-C.  Substituting  Eqs . (43a) -(43f)  in  Eq.(42),  and 
omitting  the  same  factor  as  that  in  (41)  we  have  cp  in  the  following  form, 

T’/Z  r *<1  *<0  1 

^ . J desec=>e  (I,,2"i;,2)] 


tT/2  k 

^Jdesec^e  ^(1,  ,.i;  ,) 


^'2,1'^'2,1^  " k,-k.  ^'2,2'^'2,2^] 


^A^-=U^C3.r';.,)-^C3.2S.2)] 


In  a similar  manner  we  have  in  the  following  form; 


n/2 

kj 

k 

^''l.2'^‘'i.2^] 

«Px"  ” J'  d0sec*0 
0 

Lk,-k2 

^'^l,l'^'^l,l^ 

2 

k 

K, 

: tt/2 

k, 

+ r d0sec^0 

1 

1 

^'‘2,l‘^'*2,l^ 

4 

•^3 -*^4 

^'^2  2'^'^2  2^^ 

j ec 

j + J d05ec*0 

1 0 

1 

1 

Lk3-k4 

^''3.1'"''3,P 

J-  ^ 

f 

1 

^'^3.2'''*3.2M 

Substituting  (I.  .+  ij  .)  and  (J.  .+  j!  .)  (which  are  given  in  Appendices 
' tj  • *J  • ij  • tj 

A,  B,  C)  in  the  above  equations  (44)  and  (45)  , making  use  of  Equation  (40)  and 
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omitting  the  static  elevation  and  time  factor,  the  elevation  C is 
determined  in  non-dimensional  form, 


pgC  . _ } T 

n 


2 (Uk,cose-o)^ 

^ k (k  -k  I ^ 2 

i-l  1 2^ 


V , + (Uk  .cose-t-c)^  1 


^ r'  {t^  (r:-iin 

+ sinB  ^k^k,  '1  1 ' 

> j ^ 


where  * indicates  the  complex  conjugate  , and 


^EC-ir  (logjsj  . 
m*l 


m*  I 


nvfl  '*^jt^m 


i = 1,2 


R-2040 


m=l 


-i  k,s 


V , ..m+l  3 mf-  r .,  n i2rr\  ^ , 


l£.(k,)-i;(k  ) = 

5^538  -ik.s  , 

E (-1)  e "’|E^[-ikjsJ-  oj  , 0,  Rekj<0 

m=l 


\ 


m=J 


g 

L (-1)"^’ e"' Vm{Ej[-ik^s^]+  q},  ^ 0,  Rekj^<  0 (52) 


s as  well  as  the  derivations  of  the  above  formulas  are  given  in 
m 

Appendices  A,  B,  C and  D. 


The  functions  ( I j “ I j)  i ( I 2“ ' 2^  ’ (^7)  to  (52)  contain  sums  of 

logarithmic  functions.  It  might  seem  that  each  function  has  singularities 
at  points  (a,b) , (-a,b) , (a,-b)  and  (-a,-b)  where  a pair  of  s^  vanish  (see 
Eq.A-5).  But  one  can  verify  that  the  function  has  no  singularities  when 
all  the  terms  are  taken  together  or  as  one  might  also  put  it,  any  singu- 
larities in  the  individual  terms  are  cancelled  due  to  the  alternating 
series.  Thus,  the  solution  is  regular  everywhere. 

As  a special  case,  consider  the  wave  elevation  due  to  a uniformly 
moving  pressure  distribution  without  oscillation.  Then  according  to  Tqs. 
(25b) , (26b)  and  (26c)  we  have, 

T = 0 , kpk^=0 

k2=k4='^o^®^^®  (53) 
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Hence,  Eq.(**6)  is  reduced  to 


pgC 

n 


(53a) 


and  further  since  1 1 2(1^2^  “ ’ ^ I ’4 (*^{4)  “ * 4(*^l4)  1 *''®  complex  con- 

jugates of  each  other,  the  above  is  expressed  In  the  form 


^■--’7  J • (S’") 

“ 2ti  o 

The  resultant  elevation  Is  the  sum  of  the  dynamic  and  static  terms  as 
shown  in  Eq.(40).  The  static  pressure  H is  discontinuous  at  Ixl^a  . 

The  value  of  the  discontinuity  is  defined  as 

n(ia,y)  - ^[n(±a  + 0,y)  + n(ia  - 0,y)J  . (54) 

Hence,  the  total  non-dimensional  elevation  is  obtained  by  adding  the  follow- 
ing to  the  dynamic  term  (46)  or  (53b), 

¥ -1'^ 

static  I-  1 IxKa  (54a) 


20 


R- 20^*0 


ESCAPE  AREA  AND  VOLUME 


The  escape  area  A at  the  stern  section  and  the  volume  V under 
the  plenum  area  are  obtained  by  integrating  the  wave  elevation  Q of  Eq.(ti6) 

b 

A - J Cdy 
-b 


a b 

V • J'  J Cdydx 

-a  -b 


(55) 


Referring  to  Eq.(i46),  it  is  seen  that  the  above  integrations  have  to 
use  the  following  component  integrals; 

For  the  Escape  Area  — 


”TTra('3*'3) 


-b 

-b  * 

"-Tk^sine  {■^'l'*''! 


b 


(56) 


For  the  Volume  — 

I J(l  .|j)dyd»  - - 

-a  -b  -b,-a 


a b 


I 


-a  -b 

a b 


k*  sinOcosG 


{k^(l3+lp+(l  l+l  |)-^l2(k^)  + l2(kjj)jj- 


b,a 


b,  -a 


where 


(5  7) 


<»  ki 


. Jdk  V 


(5® 
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These  are  evaluated  in  Eq.(A-33)  in  Appendix  A. 


Using  Eqs.(46),  (56)  and  (57)  t and  omitting  the  time  factor  e^'^^  , 
the  non-dimensional  escape  area  and  volume  are  written  in  the  following 
forms : 


o* 


n/2 

Jde 


t 


■i"  “4*  (nr  IV'JV  ^ 

o sin  0 ^'^3  4 ^ £-3  kjj  (kj-k^)  ^ £ 5 £JJj 


nr 


I: 


+ [de 


7k  j 

sec^0 

1-4- 

sin^0 

Ik^ 

sec^0 

sin^0 

lk3k^ 

sec^0 

sin^0 

lk3kj^ 

1 

4 

sec  0 

sin^0 

lk|k2 

sec^0 

(_£L( 

5in^0 

4 

sec  0 

sin^0 

r 2 , , (Uk,cos0-O)®  P 


£■3  kjj  (kj-k/j) 


(59) 


f^2  4q  f a 2 , . (Uk .cos0-o)^i 

J-d6i^{j2p  (I  +i')-s  i— — [k  (I  +|')+(|  ti  ) 

o sin  0 *‘r2  ° ® £«=|  ‘^£  L z J 3 II 


-(l2(k£)+i2(k£))]} 


£-3  (kj-k^j) 


[k,(.;..f).(,>,;’^). 

-(l4(k£)+i;(k£))]} 

f)+{i>i;’V 

-(l5(k£)  + i5(k£))]}| 

(60) 

In  the  above  Eqs.(59)  and  (60),  the  symbols  for  the  Integral  limits  have  been 
omitted. 


£■*3  kjj  (kj“kj^) 
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Since  pairs  of  (Eqs.A-5)  vanish  at  the  corner  points  such  as 

Sj«s^*0  at  x«a,  y*b,  etc.,  the  logarithmic  functions  and  the  exponential 
integrals  in  Ij,  ||,  l2(^£)*  dc.,  become  evidently  singular  at  the  points. 

Examination  shows,  however,  that  the  sums  i -( I ,+l  !)+l, (k*)+l.l(k  .)  [ , 

r*!*  I’l 

|-(lj  + |j  )'*‘l4(k£)  + l^(kp|,  etc.,  are  all  bounded  at  these  corner  points 
and  each  1^,  1^,  1^,  is  bounded  as  well  at  these  corner  points. 

The  regularity  of  the  above  expressions  in  Eqs.(59)  and  (6o)  are  the  ex- 
pected results  which  derive  from  the  integration  of  the  regular  function  C 
In  Eq.(46) . 

For  the  special  case  o-o'^  the  escape  area  at  the  stern  and  the  volume 
under  the  plenum  area  are  obtained  from  Eqs.(59)  and  (60)  by  using  the 
conditions  (53): 


*^^0-0  2^Lk^  o sin®0 


(6la) 


o-*0  o 


[l2(k^sec^e)  + l2(k^sec=0)]} 


-b,-a 


(61b) 

The  resultant  escape  area  and  the  volume  are  obtained  by  adding  the  static 
terms  to  the  dynamic  terms,  (59)  to  (6lb). 

Using  Eq.(54),  the  static  escape  area  at  the  stern  and  the  static 
volume  are  given  in  r.on-dimensional  form  by 

0,  - ~£aA b 


St  jjlS 


(61c) 


static 


where  a and  b are  half-length  and  beam  of  the  rectangular  pressure  dis- 
tribution and  Lisa  nominal  length  of  the  SES  or  ACV. 

’''See  foonote,  p.24. 
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The  dsymptotic  expression  of  0/  end  v for  the  frequency  o~»  are 
obtained  from  Eq,(59)  and  Eq.(60),  respectively.  Since  0-*  n/2  for  o-*» 
in  Eq.(29),  the  second  integrals  in  Eqs.(59)  and  (60)  vanish  for  a-»«  . 
Also,  the  terms  contained  in  the  summation  symbols  2 should  vanish  for 
o -CO  due  to  the  express  ions  for  the  wave  numbers  as  shown  in  Eq.(25).  The 
terms  with  * in  Eq.(59)  and  Eq.(60)  are  the  complex  conjugates  of  those 
which  are  independent  of  the  frequency  o and  the  speed  U . 

In  view  of  the  above  remarks,  we  have  the  asymptotic  expression  of 
the  resultant  cr  and  v for  o-c* 


sin*6 


In’(l3+l3)  - I 


- 0 


(62a) 


v 


r 


l__ 

2n*L* 


TT/2 

J d0 


sec^0 

810*0 


Rede+li) 


V 


i 


0 


(62b) 


where  the  suffixes  r and  I Indicate  the  real  and  imaginary  components, 
respectively.  In  the  above  integrations,  the  symbols  for  the  integral 
limits  have  been  omitted. 

From  the  expressions  of  of  and  v , it  is  seen  that  they  are  In  phase 
with  the  pressure  11  at  the  frequencies  c— 0 and  a— oe  . 


It  should  be  noted  that  in  the  limiting  cases  when  a — 0 or  o — 00,  then 
a and  t can  be  interchanged.  * 
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THE  RELATIONSHIP  BETWEEN  BUBBLE  PRESSURE  AND  HEAVE 


THe  air  pressure  in  the  plenum  of  the  SES  (surface  effect  ship)  is 

perturbed  by  the  small  vertical  oscillatory  motion  heave  ze*'^^  . The  first 

order  perturbed  pressure  p.e*^*  which  is  dynamic  and  known  as  bubble  pres- 

w lot 

sure  will  be  obtained  by  relating  it  to  the  first  order  heave  ze  in 

the  following  analysis. 

Since  the  volume  and  pressure  changes  are  assumed  to  occur  rapidly,  the 
state  of  air  pressure  and  volume  in  the  plenum  is  assumed  to  satisfy  the 
adiabatic  law, 

pv^- const.  (63) 

where  Y Is  the  ratio  of  specific  heats  which  is  1.4  for  air.  The  differ- 
entiation of  the  above  with  respect  to  time  t gives 


Expanding  p and  v in  the  foregoing  up  to  the  first  order,  we 

have 

P - p‘°>  ♦ « P<'>  ,45, 

, - ♦ . v"' 

where  p^^^  and  v^^^  are  the  pressure  and  volume  in  the  equilibrium  stage 
p^'^  and  v^*^  are  the  first  order  perturbed  pressure  and  volume,  respec- 
tively, and  e is  a parameter. 

Substituting  Eq. (6S)  into  (64)  and  collecting  the  first  terms  only, 
we  then  have 

dp<')  .p^O)  evf) 


where 


(I)  lot 

> - Pb® 


the  perturbed  first  order  volume 


■ p^  + p^  ■ total  pressure 


/o)  - 
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>■  air  pressure  in  equilibrium  or  ambient  pressure 
p • atmospheric  pressure 

O 

■ the  volume  in  equiiibrium 

Thus,  Eq.(66)  is  written  with  omission  of  the  time  factor  in  the  foi lowing 
form 

‘^b  • “"''b 

o 

where  iav^^  is  the  variation  of  the  volume  flux  which  is  positive  in  the 
direction  outward  normai  to  the  plenum  chamber  as  indicated  in  Figure  Ic 
and  is  the  sum  of  the  flux  normal  across  air-water  interface,  pienum  ceiling, 
fans  and  seal  apertures.  The  above  relation  and  the  normai  voiume  fiux  have 
been  used  in  References  3,  II,  12  which  confirm  that  the  use  of  (68)  is 
identical  to  that  of  the  method  of  momentum  change  in  Reference  I.  The  normal 
volume  flux*^  has  been  expressed  by  the  surface  integral  of  the  normai  velocity 
which  is  derived  by  the  Green's  function  technique. 

I . Over  the  Air-Water  Interface 

Equations  (60) , (6lb)  and  (6lc)  give  the  volume  due  to  the  water 
wave  elevation  per  unit  pressure,  distributed  uniformly  over  the  air-water 
interface.  Hence,  the  contribution  to  the  volume  fiux  duo  to  the  bubble 
pressure  is  given  by 


iov.  ■ -lav  p. 

b p g '^b 


(69) 


where  v is  the  non-dimensional  volume  due  to  wave  elevation  per  unit  uniform 
pressure  over  the  rectangular  patch  and  is  the  water  density.  The  minus 

sign  indicates  the  sense  of  the  flow  direction. 


2.  Across  the  Plenum  Cellinq 

The  variation  of  the  volume  flux  normal  outward  to  the  plenum  chamber 
across  the  ceiling,  due  to  its  heaving  displacement  z , Is  given  by 

iovjj  - +loA^z  (70) 

where  A.  is  the  cushion  area  or  the  area  of  ceiling.  The  -f  sign  indicates 
c 

the  sense  of  the  flow  direction. 
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3.  Across  the  Fan  Aperture 


The  variation  of  the  volume  flux  which  flows  Into  the  plenum  across 
the  fan  aperture  is 
SQf 

" BiT  Pb  ^71) 

sdf 

where  — Is  the  discharge  by  fan  per  unit  pressure  per  pair  of  fans.  The 
minus  sign  in  Eq.(71)  indicates  the  sense  of  the  flow  direction  and  n is 
the  number  of  pair  of  fans. 

4.  Across  the  Stern-Seal  Aperture 

The  variation  of  the  volume  flux  flowing  outward  from  the  plenum 
across  the  seal  aperture  in  the  plenum  wall  is  estimated  by  applying  the  per-!- 
turbation  procedure  to  the  nozzle  discharge  equation. 


k A 


p,r 


I ^2 


(72) 


* 


where  k is  the  flow  coefficient  which  varies  between  0.6  and  0.8.'  Hydro- 
dynamic  studies  have  shown  that  for  an  ideal  fluid  passing  through  a narrow 
slit  or  a circular  opening  k~  0.611  which  has  been  experimentally  verified. 


A is  the  area  of  the  nozzle  throat  or  seal  opening 
Pj  is  the  pressure  in  the  upstream  or  in  the  plenum 

P2  is  the  pressure  at  the  nozzle  or  seal  and  identical  to  the  atmospheric 

pressure  p 

3 

p_  is  the  air  density 
d is  the  discharge  . 

To  evaluate  the  increase  of  volume  flux  from  its  equilibrium  condition,  we 
expand  the  terms  di  A and  P^-P2  up  to  the  first  order  perturbed  terms  assum- 
ing that  the  terms  are  functions  of  e . Namely,  we  have  the  following: 

Q ■ + e q(*) 

A,  - A + e A^’^ 

s o 

P1-P2  " Pq  ® P^*^  (73) 

where  d^.  A^  and  p^  are  the  zeroth  order  term,  whereas  d^'^*  A^'^  and  p^'^ 

'’''As  the  "nozzle"  is  formed  by  the  water  surface  and  the  seal,  the  coefficient 
k will  depend  on  Froude  number  and  on  Reynolds  number  for  small  opening 
areas,  A , as  measured  at  zero  speed 


i 
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are  the  first  order  terms.  The  first  order  terms  are  redesignated  by 

i i pO-v'"  • 

Substituting  (73)  into  (72)  and  making  use  of  the  following  expansion 
+ 


(75) 


we  have  the  expression  for  Q,  from  which  the  first  order  equation  is  obtained 
in  the  following  form 


ik  * *») 


(76) 


The  variation  of  the  seal  opening  area  A in  (76)  is  due  to  heave  z 


and  bubble  pressure  p^  . Namely, 
6A. 


3r 


z + 


Pb 


(77) 


where 


dAs 

^ — ■ the  rate  of  change  of  seal  opening  due  to  seal  elasticity 
per  unit  heave 

dAj 

Y— ■ the  rate  of  change  of  seal  opening  due  to  seal  elasticity 
P per  unit  pressure 

-o  ■ the  rate  of  change  of  opening  due  to  the  wave  elevation  per 
Pw^  unit  pressure,  In  way  of  the  stern  seal.  a is  the  non- 

dimensional  escape  area  which  has  been  given  in  Eqs. (59) , (6la) 
(6ic),  and  (62a). 

Substituting  (77)  into  (76),  we  have  the  following  volume  flux  which  flows 
outward  in  the  normal  direction  to  the  plenum  wall  across  the  seal  aperture: 


lav. 


|2p..  r /a  6A  , \ bA,  1 


(78) 


Thus  the  itemization  of  the  variation  of  volume  fluxes  has  been  completed. 


28 


R-20i*0 


The  sum  of  the  first  order  volume  fluxes,  Eqs.(69),  (70),  (71)  and  (78), 
Is  now  substituted  Into  Eq.(68).  Solving  for  the  unknown  bubble  pressure  p^, 
we  have  the  non-dlmens lonal  expression  in  the  following  form 


Indicates  the  bubble  (or  cushion)  height. 

dA  dA 

s s 

The  dynamic  property  of  the  stern  seal  and  is  unknown.  At 

the  present  stage  of  Investigation,  we  assume  the  seal  Is  rigid.  This 
assumption  will  evidently  simplify  the  aero-hydrodynamic  analysis.  In 
addition  to  this,  the  analysis  will  be  practically  applied  to  the  model  test 
which  uses  a semi-rigid  seal,  e.g.,  in  the  case  of  the  XR-5  SES  model.  Refer- 
ence 3 reports  the  investigation  of  the  seal  stiffness  on  the  aero-hydro- 
dynamic  bubble  pressure  of  a SES  model. 

The  above  assumption  reduces  the  expression  of  the  bubble  pressure  to 
the  following  form, 


If  the  water  wave  induced  volume  v and  area  a are  set  equal  to  zero,  the 
above  relation  will  represent  the  bubble  pressure  due  to  heave  of  an  air 
cushion  vehicle  above  the  ground. 


Although  Eq.(80)  is  derived  by  means  of  perturbation  analysis  keeping 
only  the  first-order  term,  two  limiting  values  of  the  bubble  pressure  under 
extreme  conditions  can  be  obtained; 

a)  When  translation  velocity  U and  frequency  cr  approach  to  zero  with 
openings  closed  aft  and  no  fans  in  operation  (closed  box),  then 


I im  ^ 

U-0  - Pw9  ’ 

0-0  K 


(81b) 
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b)  When  the  heaving  vehicle  is  in  translatory  motion,  Uj^O,  with  all 
the  openings  in  effect,  then  when  frequency  a - 0 

-«  0 with  90°  phase  lag.  (8 lb) 

O-*0  2 

K 

where  the  following  information  has  been  used,  namely,  for  a “•  0 , accord- 
ing to  (59)  to  (61); 

1)  and  v^  are  negative  real  constant  values 

2)  (Xj  and  v.  are  linearly  proportional  to  a 

The  above  relation  is  also  equally  valid  for  the  uncoupled  aerodynamic 
analysis,  namely,  for  v=Of=0  . 


It  is  apparent  that  the  asymptotic  expression  of  the  bubble  pres- 
sure Eq.(80)  for  a OB  is  dependent  on  the  behavior  of  the  wave- induced 
volume  V for  o -•  ob  on  the  air-water  interface.  Since  the  asymptotic  expres- 
sion of  V , Eq.  (62b),  is  v^  = a negative  constant  and  v.  = 0 , we  have 
with  water  wave  effects 


(82a) 


whereas  with  no  water  wave  effects, 


lim  ^ /z  t=  _y(  1 + ^)  . (82b) 

a-oopo/  h \ pj 


Both  limiting  values  of  pressure  signatures  show  a l80-degree  phase  lag. 


*|t  should  be  noted  that  in  the  limiting  cases  when  a - 0 or  a -*  ®,  then 
a and  t can  be  interchanged. 
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ADDED  ^WSS  AND  DAMPING  COEFFICIENTS 

The  real  and  imaginary  components  of  the  bubble  pressure  Eq.(80)  are 
given  by 


p 1 ■ 

-'<(l■>Pa/Po)Ll-V(l  * P./Pp)  -^1 

' [l-V(l.p,/p„)  *[v(l*P,/p„) 


■^Ttn 


2 '’o«i 


- C''(l»Pa/b„J  -Vo 


r Po^rf 

r Po^n 

•-  o -* 

{y(Hp3/Po) 

where 


. k _ L 


Q = V = + — Of.  

>■  »■  Pw9  ^ ' Pw9  ^ Pa 

«i  ■ ’i  ^ p ("Ivl  *'2p:  ' “r 


<1-  " V.  r — 
' ' P 9 


The  integration  of  the  bubble  pressure  (83)  over  the  cushion  area 
yields  the  aero-hydrodynamic  force  due  to  heaving  motion  of  a surface 
effect  ship  in  translation  on  the  calm  water  surface.  The  real  and  imaginary 
comoonents  of  the  force  are  designated  as  aero-hydrodynamic  inertial  and  damp- 
ing forces  from  which  the  added  mass  M"  and  damping  factor  N are  defined 
in  the  following  form, 

Pbj.^c 


Pb  ^c 


where  the  minus  sign  comes  from  the  Re(pjjA^e  )• 
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The  non-dimensional  added  mass  and  damping  coefficients  of  Eqs.(83) 

7 


can  be  defined  in  the  following  form, 
M"  . . ^ N 


'aT  . 
2 Pw*- 


; c . 

3 a 


I 


(86) 


The  added  mass  M"  and  damping  factor  N for  o -*0  are  derived  by 

knowing  the  behavior  of  the  a , a. , v and  v.  at  a -•  0 ; 

r I r I 


O " (W  ^ ® 


dv. 


> 0 


j da 

'r-'IM  * 


or 


IM 


(87) 


The  results  are  given  by 
.a 


M" 

a-»0 


E A‘ 
D‘ 

.3 


N * -r^ 
a—  0 0 


(88) 


where 


a—0 


P 9 


:'o ^ 


(89) 


where  D is  always  positive,  whereas  E may  take  positive  or  negative  value. 
Thus  the  damping  factor  N is  always  positive  and  the  added  mass  may  be 
positive  or  negative  (i.e.,  the  motion  of  the  vehicle  and  the  added  water 
mass  are  in  opposite  directions).* 


Definite  signs  are  obtained  when  the  water  wave  effects  are  omitted. 
Name I y , 

(90) 


^(Pa-^Po> 
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From  Eqs.(90)  we  have  always  the  negative  added  mass  and  positive  damping 
factor  for  the  SES  without  water  wave  effects  or  ground  effect  machine 
at  a -»  C. 

It  is  apparent  from  Eqs.(88),  (89)  and  (90)  that  the  added  mass  M" 
and  damping  factor  N at  the  vanishingly  low  frequency  are  proportional 
to  the  cushion  area  squared  and  dependent  on  the  contributions  by  the  fan, 
seal  leakage,  wave- induced  area  and  volume  and  the  cushion  volume  and 
ambient  pressure  which  are  chosen  in  the  preliminary  design  stage. 
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SCALE  EFFECT  ON  THE  BUBBLE  PRESSURE 


Reference  12  presents  the  scale  effect  on  the  real  component  of  the 
bubble  pressure  In  the  two-dimensional  SES  against  the  scale  ratio  \ ■ 
ship  length/model  length.  It  shows  that  the  scale  effect  is  larger  for 
the  low  frequency  than  for  the  high  frequency. 

A similar  study  can  be  carried  out  for  the  present  three-dimensional 
SES  model . 

Since  ^o^i^'^o  dimensionless,  the  scale  ef- 

fect on  the  magnitude  of  the  pressure  (or  the  components)  can  be  estimated 
by  scaling  the  ambient  pressure  p^  by  \p^  . For  the  real  component  the 
ratio  of  the  coefficient  for  ship  and  model  is 


^ship  r 


Pk/P. 


model 


1 + Pg/^pj 

I-Y(I 

- 

+ Pa/^Po>Po‘ir^V 

[l-Yd  + Pa/^Po)PoQ/v 

'-13 

oJ 

LY{I  + Pa''^Po^Po^i''V 

a 

A'  Pa/Po)[' 

-Y(I 

+ p /p  )p  Q /v 
^3  ^ ~ o^r  o 

_1-V(I  + Pa/Po)Po'ir/V 

'd 

+ 

+ Pa/Po^Po^i^'oT 

(90 


Similar  expressions  can  be  obtained  for  the  imaginary  component  as  well 
as  for  the  resultant  magnitude. 
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NUMERICAL  CALCULATIONS  AND  DISCUSSION 
Water  Waves  - Elevation,’  Escape  Area  and  Volume 

The  t'lWdel  used  In  Reference  2 has  been  chosen  to  compare  the  analytical 
values  of  wave  elevation,  escape  area  at  stern  and  volume  on  the  air-water 
interface  in  the  developments  of  that  reference  and  the  present  study. 

A scries  of  calculations  have  been  carried  out  for  the  elevation  of 
the  free  surface  along  the  wall  y*b  at  Froude  number  F^^O.S  for  different 
values  of  frequencies  T-0,  and  t-0.2375  and  T»0.2750  which  bracket  the  criti- 
cal frequency  T"0.25.  The  results  are  shown  in  Figures  2a  and  2b  and  compare 

v'ery  favorably  with  those  of  Doctors  who  defines  Fp  ■ v/JUag  , Figure  2c  is  i 

! 

'is  a comparison  of  the  va'ues  of  wave  elevation  of  the  corner  point  x/a»-l  T 

and  y/b»l  in  terms  of  frequency  t,  evaluated  by  the  present  analysis  and  by  1 

f 

that  of  Reference  2.  In  the  low  frequency  regime  the  agreement  is  excellent,  f 

whereas  outside  that  the  discrepancy  is  quite  substantiai  especially  in  the 
imaginary  part  of  the  elevation  values. 

Similar  calculations  were  carried  out  for  the  escape  area, a , at  the 
stern  and  the  volume,  v , on  the  air-water  interface  for  the  entire  frequency 
range  of  practical  interest.  A comparison  of  the  results  of  both  analytical 
solutions  is  shown  in  Figures  3a  and  3b  for  the  escape  area  coefficient,  a , 
and  volume  coefficient,  v , respectively,  for  a series  of  frequencies  t . 

The  agreement  is  rather  satisfactory  up  tot ’•1.5  except  for  the  real  value 
of  the  V coefficient.  It  Is  interesting  to  notice  that  the  value  of  the 
escape  area  oscillates  with  t,  whereas  the  volume  coefficient  v behaves 
relatively  monoton  I ca 1 ly.  Sudden  changes  of  the  values  of  of  and  v are 
exhibited  at  the  critical  t-0.25  as  should  be  expected.  Furthernrare,  it  is 
interesting  to  note  that  as  t-»0  , the  imaginary  part  of  both  quantities  tends 
to  zero  whereas  the  real  parts  tend  to  finite  negative  values.  These  limiting 
values  can  be  readily  obtained  from  Equations  (61a),  (61b),  and  (62a),  (62b). 

As  the  frequency  approaches  higher  values,  the  imaginary  parts  of  a and  v 
diminish  showing  a trend  to  zero  asymptote,  whereas  the  real  components  tend 
to  a negative  non-zero  asymptote.  The  non-zero  values  are  less  than  the 
static  values  of  Ofj."-0.25  and  v^"-0.25.  It  is  apparent,  therefore,  that 
the  dynamic  values  of  a and  v in  Eqs.(62a)  and  (62b)  are  positive  and  their 
magnitudes  are  less  than  their  corresponding  static  values  (Eq.6lc). 
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similar  calculations  have  been  carried  out  for  XR-5  SES  model  at 
Fn*  0.72,  for  which  extensive  experimental  work  is  avai lable.^”^  Figure  4 
is  a sketch  of  the  model  and  Table  I presents  its  principal  particulars. 

The  results  for  a and  v are  exhibited  in  Figures  5a  and  5b  versus 
the  frequency  t . The  trends  of  the  results  of  both  models  are,  in 
general,  similar  except  that  the  oscillatory  behavior  of  a has  been  re- 
duced considerably  in  the  second  model. 

Bubble  Pressures 

The  values  of  bubble  pressures  in  the  plenum  calculated  from  the  ex- 
pressions given  by  Eqs.(80),  (8I)  and  (82),  for  the  XR-5  SES  model  at 
Fn"0.72  are  graphically  exhibited  in  Figures  6a  and  6b  for  a scries  of 
reduced  frequencies  t . 

Figure  6a  shows  the  values  of  bubble  pressure  and  phase  shift  angle 
versus  T in  the  presence  or  absence  of  water  free  surface  effects.  In  the 
limit  as  T-‘0  , the  expression  given  by  Eq.(8l)  shows  the  pressure  ampli- 
tude reduces  to  zero  with  phase  lag  of  90  degrees.  The  pressure  amplitude 
generally  increases  as  T increases  tending  to  the  asymptotic  value  given 
by  Eqs.(82a)  and  (82b)  for  both  cases,  with  and  without  the  free  surface 
wave  effects. 

Figure  6b  gives  the  entire  picture  of  bubble  pressure  versus  T ranging 
from  zero  to  6OO.  It  is  seen  that  the  amplitude  reaches  its  asymptote, 

V(l+  P^/Rq)  at  t«  200,  whereas  the  phase  angle  may  reach  its  asymptote  at 
t"800.  It  Is  to  be  noted,  however,  that  the  values  of  the  bubble  pressure 
in  the  uncoupled  condition  (i.e.,  without  water  wave  effects)  are,  in  general, 
much  larger  than  those  obtained  under  the  coupled  condition.  In  fact,  the 
asymptotic  value  of  the  former  is  about  thirty  times  larger  than  the  latter. 

Added  Mass  and  Damping  Coefficients 

The  heave  added  mass  and  damping  coefficients  C , C . of  the  XR-5  SES 

0 d 

model  have  been  calculated  and  are  shown  in  Figures  7^  and  7b.  The  experi- 
mental values^  for  the  model  are  also  compared  with  the  theoretical  values 
In  Figure  7a.  Both  results  show  the  same  trend  and  the  experimental  damping 
coefficients  agree  better  with  the  present  theory  when  wave  effects  are  taken 
Into  account.  The  asymptotic  values  of  the  coefficients  obtained  from  Eqa.(85) 
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and  (88)  as  T-*0  are  the  same  for  both  coupled  and  uncoupled  conditions. 

The  Imaginary  part  of  the  bubble  pressure  from  which  the  dampling  factor  N 
Is  determined  is  linearly  proportional  to  the  frequency  t (or  a)  at  low  fre- 
quency for  both  coupled  and  uncoupled  conditions  (see  Eq.85).  The  experi- 
4 

ments  also  find  this  component  linearly  proportional  to  the  frequency  for  an 
ACV  model  over  the  ground.  The  added  mass  coefficient  is,  in  general,  negative 
for  the  specific  model  except  in  a narrow  region  as  indicated  in  Figure  7a. 

The  non-dimensional  coefficients  and  C.  diminish  as  T Increases. 

d d 

As  T approaches  « , the  coefficients  tend  to  vanish  for  both  coupled  end  un- 
coupled conditions  (see  Figure  7b). 

Scale  Effect  on  the  Bubble  Pressure 

The  scale  effect  on  the  bubble  pressure  amplitude  has  been  calculated 
for  the  XR-5  SES  model  at  Fj^*  0.72,  and  the  results  are  shown  in  Figure  8. 

It  Is  seen  that  the  scale  effect  is  increasing  with  increasing  frequency  at 
constant  scale  ratio  \ and  increases  also  with  \ at  constant  frequency. 

The  scale  effect  increases  with  both  increasing  frequency  and  scale  ratio.  For 
instance,  the  ship  pressure  at  t»4.5  is  underestimated  by  I4%  when  the  model 
scale  ratio  \*10  , which  corresponds  to  about  a 133  long-ton  ship* 

12 

The  previous  work  on  the  scale  effect  is  concerned  with  the  real  com- 
ponent of  the  pressure  only.  It  shows  that  the  scale  effect  is  larger  at  low 
frequency,  a trend  which  Is  opposite  to  the  present  result.  However,  it  is 
to  be  noted  that  the  analysis  of  Reference  12  deals  with  the  two-dimensional 
theory  and  hence  it  is  difficult  to  compare  with  the  present  three-dimensional 
analysis. 


Heave  Resi 


Since  the  theoretical  heave  added  mass  and  damping  coefficients  are 

available,  the  heave  response  with  pitch  restrained  has  been  calculated  by 

6 

utilizing  test  results  for  the  restoring  and  exciting  forces  of  the  XR-5  SES 
mode  I . 


Figure  9a  shows  the  mean  line  drawn  among  the  experimental  values  of 
the  non-dimensional  heave  exciting  force^  versus  the  non-dimensional  encounter 
frequency  cVtTg  . The  extrapolation  has  been  made  by  assuming  that  the 
heave-exciting  is  identical  to  the  restoring  force  at  zero  frequency.  The 


*These  conclusions  reached  after  the  set  of  calculations  performed  for  the 
XR-5  SES  for  given  fan  characteristics  and  discharge  should  be  used  with 
caution  and  cannot  be  generalized  before  extensive  calculations  are  performed. 
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mean  line  has  been  cut-off  at  the  frequency  beyond  which  the  data  are  un- 
reliable. 

The  experimental  added  mass  and  damping  coefficients  have  been  obtained 
from  the  three  test  points  indicated  in  Figures  I8  and  19  of  Reference  6. 

Using  these  values,  three  heave  responses  have  been  calculated.  Figure  9b 
compares  the  semi-analytical  and  the  experiment-based  calculations.  The  semi- 
analytical  response  is  higher  than  the  experiment-based  calculation  results 
due  to  the  fact  that  the  theoretical  damping  coefficients  are  lower  than  the 
experimental  values. 

The  theoretical  calculations  show  a similar  trend  to  that  of  the  heave 

with  pitch  free  (x  symbols)  which  is  obtained  from  Reference  7.  The  natural 

pitching  frequency  of  the  model  is  5 rad/sec.  It  is  seen  in  Figure  9b  that 

the  theoretical  heave  at  the  natural  pitching  frequency  is  about  50%  of  the 

heave  with  pitch  free.  This  kind  of  characteristic  heave  response  in  coupled 

8 9 

and  uncoupled  conditions  for  an  SES  model  has  been  found  by  Fridsma  ’ at 
Davidson  Laboratory. 

Numerical  Procedures 

The  computational  work  has  been  carried  out  at  the  Stevens  Institute  Com- 
puter Center  on  the  POP-IO  high-speed  digital  computer.  Program  structure 
with  Flow-Chart  of  the  main  subroutine  VLLI,  program  listing  and  descriptions 
of  input  and  output  are  given  in  Appendix  E. 

The  exponential  integral  of  complex  argument,  which  is  the  basic  element 
in  all  the  analytical  expressions,  is  discussed  in  detail  in  Appendix  0 with 
the  view  to  utilizing  high-speed  digital  computers. 

Furthermore  for  convenience  in  the  numerical  integration  of  Eqs.(46) 

and  (59)  to  (62),  a transformation  is  introduced  to  simplify  the  structure 

of  the  singular  behavior  at  0 (critical  point,  see  Eq.29),  and  then  the 

c 

"Lagrange  interpolation  method"  has  been  employed  for  the  evaluation  of  the 
finite  contribution  of  the  high-order  singularity. 
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SUMMARY  AND  CONCLUSIONS 

A mathematical  model  has  been  developed  for  the  evaluation  of  the 
added  mass  and  damping  coefficients  of  a heaving  surface  effect  ship  !n 
uniform  translation  on  a calm  water  surface. 

The  analysis  furnishes  the  variation  of  bubble  pressure  in  the  plenum 
and  its  scale  effect  in  terms  of  reduced  frequency  t in  the  entire  range 
of  practical  interest,  and  also  provides  information  on  the  wave  elevation 
of  the  free  surface,  the  escape  area  at  the  stern  and  the  volume  on  the  air- 
water  interface. 

A program  has  been  written  in  FORTRAN  IV  language  and  applied  to  the 
PDP-10  high-speed  digital  computer. 

The  following  basic  assumptions  were  made: 

1.  The  deformation  of  the  water  surface  participates  in  the  generation 
of  the  bubble  pressure  in  conjunction  with  the  action  of  seals  and  fans. 

2.  The  deformation  of  the  free  surface  under  the  oscillatory  rectangu- 
lar pressure  patch  in  uniform  translation  is  used  to  display  the  way  in  which 
the  motion  of  the  water  surface  participates  in  the  evaluation  of  the  pressure 
variations  in  the  plenum  air. 

3.  In  the  dynamic  analysis  of  the  air  flow  phenomena,  pressure  and 
volume  changes  are  assumed  to  occur  sufficiently  rapdily  so  that  the  adiabatic 
law  governs  the  basic  thermodynamic  variations  in  the  air  cushion. 

4.  The  bubble  (or  cushion)  pressure  is  assumed  to  be  spatially  homo- 
geneous. 

5.  The  leakage  of  the  air  takes  place  only  at  fan  and  stern  apertures. 

6.  The  perturbed  (or  bubble)  pressure  and  air  volume  have  the  same  order 
of  magnitudes  as  the  heave  displacement. 

Expressions  for  the  wave  elevation,  escape  area  and  volume  have  been 
derived  and  their  limiting  values  at  T-»0  and  T-*®  have  been  studied. 

The  linear  equations  of  bubble  pressure  and  heave  motion  have  been  de- 
veloped, applying  perturbation  procedure  to  the  adiabatic  law,  and  their 
asymptotic  values  at  t-*0  and  have  been  established. 
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A series  of  numerical  calculations  for  the  SES  model  have  been  carried 
out  for  the  wave  elevation,  escape  area,  volume,  bubble  pressure,  added  mass, 
damping  factor,  the  scale  effect  on  the  bubble  pressure  and  the  heave  response 
with  pitch  restrained. 

All  the  above  results  have  been  compared  with  available  experimental 
data  and  their  behavior  as  functions  of  the  frequency  discussed. 

The  following  conclusions  can  be  drawn  from  this  investigation; 

1.  The  assumption  of  a rigid-free  surface  is  not  valid.  Considerable 
differences  have  been  observed  with  and  without  the  free  surface  waves. 

2.  The  bubble  pressure  in  the  plenum  has  zero  amplitude  with  90^  phase 
lag  as  T-'O,  whereas  at  t-*“  it  reaches  a constant  value  with  |80°  phase  lag. 

3.  The  added  mass  and  damping  factors  as  T-*0  reach  asymptotically 
finite  values  whereas  as  t-*®  they  vanish. 

4.  The  scale  effect  on  the  bubble  pressure  increases  with  increasing 

frequency  and  scale  ratio.  (More  extensive  calculations  are  needed  before 
generalizing  this.) 

5.  The  present  method  can  be  applied  to  predict  the  heave  response  with 
restrained  pitch  motion  utilizing  experimental  values  of  the  restoring  and 
exciting  forces. 

6.  The  added  mass  and  damping  factor  compare  favorably  with  the  experi- 
mental data.  Unfortunately  a very  limited  number  of  experiments  are  at  present 
available  and  before  more  experiments  are  conducted,  no  general  conclusion  as 
to  the  accuracy  of  the  theoretical  prediction  can  be  drawn. 

It  is  expected  that  the  present  method  can  be  extended  to  calculate 
I)  the  heave  exciting  force  of  a captive  SES  model  uniformly  running  in  head 
regular  waves  and  2)  the  heave  restoring  force  acting  on  an  SES  model  in 
uniform  translation. 
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TABLE  I 

PRINCIPAL  CHARACTERISTICS  OF  XR-5  SES  MODEL 


Length  overall,  ft 

Length  of  bubble,  ft 

Beam  of  bubble,  ft 

Height  of  bubble,  ft 

Test  displacement,  lb 

Nozzle  factor 

Fan  flow  rate  per 

pair  of  fans,  ft® /lb-sec 

No.  of  pairs  of  fans 

Leak  area  at  equilibrium 


Froude  number 


Encounter  frequency,  rad/sec 
Natural  pitching  frequency,  rad/sec 
Heave  restoring  force,  Ib/ft 


LOA  - 15.58 
L - 13.83 


b - 2.12 


0.03 


A^°^-  0.333 

Fn-  0.72 

lu  • 3~  22 
e 


B = 1350 


INAL  SECTION  y/b= 


FIGURE  2b.  ELEVATION  ALONG  THE  LONGITUDINAL  SECTION  y/b 
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FIGURE  3a.  ESCAPE  AREA  COEFFICIENT  or  AT  STERN  FOR  THE  MODEL  IN  REF. [2] 
VERSUS  T FOR  F„  = 0.5 


FIGURE  3b.  VOLUME  COEFFICIENT  v VERSUS  T 

FOR  Fn  = 0.5  FOR  THE  MODEL  IN  REF.  [2] 


AND  MIDSTATION  SIDEWALL  SECTION 


FIGURE  5a.  ESCAPE  AREA  COEFFICIENT  a AT  STERN  VERSUS 
AT  F_  = 0.72  for  XR-5  SES  MODEL 


AT  = 0.72  FOR  XR-5  SES  MODEL 


Phase  Lag.  deg 


FIGURE  6a.  BUBBLE  PRESSURE  COEFFICIENTS  VERSUS 
AT  F„  = 0.72  FOR  XR-5  SES  MODEL 


Asymptotic  Phase 
Due  to  Eqs.(82a)  & (82b 


FIGURE  6b.  BUBBLE  PRESSURE  IN  THE  HIGH  FREQUENCY  REGION 


A Experimental  Added  Mass  & Damping  Coefficients 
— No  Water  Surface  Effects 


FIGURE  7a.  ADDED  MASS  AND  DAMPING  COEFFICIENTS  VERSUS 
AT  Fn  = 0.72  FOR  XR-5  SES  MODEL 


With  No  Water  Wave  Effects 


FIGURE  8.  SCALE  EFFECT  ON  THE  BUBBLE  PRESSURE  VERSUS  SCALE  RATIO 
AT  F » 0.72  FOR  XR-5  SES  MODEL 


1 


V 


R>2040 


Ref .[6 j 


FIGURE  9e-  HEAVE-EXCITING  FORCE  OF  XR-5  SES  MODEL  AT  F„»0.72 


* 


I 


FIGURE  9b.  HEAVE  TO  WAVE  AMPLITUDE  RATIO  OF  XR-5  SES  MODEL 
AT  Fn  - 0.72 
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APPENDIX  A 


CALCULATION  OF  I , , . I ! . , J , . . J ' (j=l,2) 

' ‘J  ' •J  ' >J  I .J 


Integration  of  I,  . and  I I in  Eqs.(43a)and  (43b)  with  respect  to 
' » J ' > J 

§ and  H yields 


I S V - - 

'l.j  sine 


[ dk 


kzi  ikSj  iks^  iks,  **^^4') 
e Ne  -e  +e  ^ ~e  • 


k(k^ 

J 


o ji  e 

^ sinDcosO  ^ 
o 


k (k-k^.) 


(A-l) 


I 


l,j  sine  ^ 


kz,  iks_  ikS/  iks^  'kSo, 

e \-e  ^ -.e  ^ -e  ^ 

k(k-kj) 


_o eK. 

sinOcosB  J 


k^(k-k.) 


(A-2) 


Differentiation  of  I,  . (A-l)  and  l!  . (A-2)  with  respect  to  x yiel( 

^ » J ^ • J 


J . = J dk 

1,1  I s I n0  ‘I 

'J  n 


^ kz.  iks,  iks,  iks  iks.N 

U cos9  r ji  ® 'e  -e  +e  ^ -e 


- T — I dk  

I s I ne  J 
o 


k{k-k.) 


.*  .=  -iL£2i|P  dk 

•j  ' stne  J 


kzi  iks-  iks,  iks.^  '*^^8) 
; N-e  +e  -e  ' +e  r 
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k(k-kj) 
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whore  J 


Si.  ■ 


Sf  - 


I , 2 iirul 

(x  cost)  + y sInO)  - (a  losO  + b sinO) 
(x  cosO  y sinO)  - (a  r.osO  - b sInO) 
(x  cosO  + y sin0)  + (a  i,osO  + b sinO) 
(x  COS0  + y sinO)  + (a  cosO  - b sinO) 
(x  cosO  - y sinO)  - (a  cosO  + b sinO) 
(x  cosO  - y sinO)  - (a  cjsO  - b sinO) 
(x  cosO  - y sinO)  + (a  cosO  + b sinO) 
(x  cosO  - y sinO)  + (u  cos0  - b sinO) 


(A-5) 


S i nee 


1 I 1 
— + 


rTkikjT  kj  k kj  k-kj 


1 


I I 


' 1.  ' ' 


Equations  (A-1)  through  (A-b)  arc  expressed  in  ttu*  followiny  forms, 
U 


I ^ I o 

'•j  VjSinO  k j'^si  nOcosO^  ''I  '2'  kjSinOcosO  '3 


^'r'2^  ■ k.sinOci 


l.J 
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kjSinO  kj*'sinGcosO' 


kjSinOcosO  ^ 


> 'isinO  ikjSinO'  2 


' . J I k j 5 i nO 
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I - j dk 


i;  - J dk 

' O 


i;  - J 


!•  . f dk 
^ o 


kz/  ikSj  iks^  'ks,  '*'^4') 
:•  'C  -e  +1;  ^ -e 


kz/  ikSr  iks,  iks,  ikSt 
. Ue  ^ .^e  ^ -e  ^ 
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Integration  of  | ^ , |^,  | J , 


To  transform  kz  + ikSj,  kz  + iks2,  etc.,  we  write 
C * X COS0  + y s i n9  - i z 

I 

Q “X  cose  - y sine  -iz 
Cj  = a cose  + b sine 
^2  “ 3 cose  - b sine 


and  we  have 


z + 

II 

I 

z + 

is5  «=  i(C'-C,) 

z + 
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Substituting  the  above  in  Ij,  1^,  ij,  we  have 
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• /riN  T r 
i,i‘  > ■ J tL'" 
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-JfL' 
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(A-14) 

1 (A-15) 


consider  the  first  integral  1,(0.  1,(0  in  Eq.(A-12)  is  analytic 

in  the  lower  half  plane  and  1,(0- 0 for  i - . Hence,  differentiating 
with  respect  to  C . we  have 


dl 


_ . i j dk 


ik(c-c,)  _^ik(;.t,)  ^^■.k(W,) 


Integration  with  respect  to  k is  given  by 


dl 


1 


1 ( > _L_  4.  J_  - 1 

■ t-Cj  ^■^'^2  ’ 


(A-17) 


(A- 18) 


Integrating  with  respect  to  C • we  have 
I ^"^^1  \ 

'l  “ 

where  K “ integration  constant. 

Since  both  1,  in  (A-12)  and  the  logarithmic  terms  in  (A-l8)  vanish  for 
^-•ioo  , we  have  K^O  . 

Hence, 


1 


C"C,  C+£,' 


x^rJ 


(A-19a) 


and 


l,|  = -(log|s,|+  i°n)-Oo9|«3l+  '°9| 

' - (A- 19b) 


z-0 


A similar  procedure  is  applied  to  1^,  I,  and  1^; 

C+C,  C-C2 

1 3 ’ '1  ^1  ^2  C+;2^ 


a4 


(A-20a) 
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«3|  --iC-(x  cose  + y sine)lj  +C,  | ( log  | | + ?^)  •.(^og|  s , 1 + p^)  1 

2^0  2“0 

+C2i(loglS2|  + iO)-(log|Sj^|  + iO^)  l3  . 


s < 0 

tn 


(A-20b) 


(A-2>a) 


‘1  “^'°^h5l'"?TT)-"('°9|"7l‘"°r.)-(’°glS6rin)-(’°9|S8|  + P„),s^^  0 (A-2»b) 


, • , C'-C,  C'+C,^ 

I3  “ I,  + £,iog  ^77^  + Tnqj 


I3  --i[ix  cose  -y  sine)  I I +C,  1 ( log | s^ | + .0) -( log | s^|  + 0^)  1 
z*«0  2=0 


(A-22a) 


+;2lOog|S8|+0^)-(log|s^|  + .0^)l], 


Integration  of  I2  and  1^ 


s^  ^ 0 (A-22b) 

m ^ 


*2  (A-8)  and  (A-9)are  identical  with  each  other  in  form;  therefore  it  is 
necessary  to  consider  the  integral 


» ik(s,-iz) 

= -J' 

O J 

_fV. 


(A-23) 


Transforming  the  above  by  translation 

w - k-k. 

J 


(A-24) 


and  rotation 


u ■-iw(5,-i2) 


(A-25) 


we  have 


A5 
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- i z) 

ik.s,+k,z  ' -u 

G = ■>  ’ J J Y 

ikj  Sj+kjZ 

The  above  is  expressed  by  the  exponential  integral, 
ik,s,+k,2 


17 


(A-26) 


(A-27) 


Hence  for  z=0, 


<2  ' " i2n  } ^ * (^-28) 


{E,[ikjsJ 


- ° 1 
i2'n  J 


for  s < 0 
m 


where  j = 1 , 2 . 


CALCULATION  OF  1^  AND  1^ 


(A-29) 


We  calculate  the  following  integrals  which  are  used  for  evaluation 
of  the  volume  in  Eq.(59): 


= J ^ 


o k 


1^  . J ^kz(-e"'''  «'''*»  -e"'"’  ) 


(A-30) 


o k 


The  same  procedure  as  that  used  previously  in  the  Integration  of  I j , 
I J , 1^  and  I J is  used. 

The  1^  and  1^  are  similarly  expressed  in  the  following  form; 

, (;)  - .eiMC+Ca)^ 

o k^ 


.,ik(C'Hx) 


o k 

For  C'“*  1“  t *6  -•  0 , and  analytic  in  C and 
I ' 

® A6 


A 


r 
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Differentiating  lg(G)  and  l^(C')  three  times  with  respect  to 
C and  G'  I and  integrating  them  with  respect  to  k , we  have 


d^, 


6 J 1 , i 1 

' G-G,  ■ G-G,  G+G,  ■ G^ 


dG 


'6 


•1 


dG 


.3 


1 . i 1.1 


(A-32) 


Integration  of  the  above  derivatives  with  respect  to  Q and  G'  yields  each 
time  a pair  of  integration  constants.  The  first  two  pairs  of  constants 
vanish,  whereas  the  last  pair  results  are  i(Gj^-G2^)  ®'^d  -iCGj^-Gj^)* 
Hence,  the  pure  constant  terms  in  lg(G)  and  Ig(G')  disappear  and  the 
and  are  finally  expressed  in  the  following  form; 

'6j  ■ ’ 7 (’‘cose+ysine)^  Ijj  -(xcos0+ysine)G,[(log|Sj  l + .O^^-'^logls^l+O^j^ 
z-0  z“0 

+{ xcose+y s i n0)  C2l_(  1 og | S2 1 +pj“(  ’ °9 1*4 1 n ) ] 

+ I [t|HO°g|*il  + t°TT)+  (’°9|*3l'"in)}  ■ 

-G2"{(log|S2|  + pJ+(log|Si^l+pJ}]  for  s^  > 0 

'sj  “ - |7(>'cos0-ysin0)^  lj|  -(xcos0-ysin0)G,[(log|S5|+9j-(log|s^j  + .o)] 

+(xcos0-ysin0)G2[(loglSg|+oJ-(log|S8|  + .o)] 


z-0 


z-0 


*6l*iiiX'‘>9| 'srin)}]  I 


for  s ^ - 
m < 0 


(A-33) 
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APPENDIX  B 


CALCULATION  OF  1^  . ;J»3,4) 

Integrating  1^  ^ and  1^  jj  with  respect  to  § and  T1  in  (A3c)  and 
(l<3d)  , and  making  use  of  (A-5)  and  (A-6)  , we  have 


, u 

o 

/ , 1 V ^ 

0 

2.-^  “ 

( k.sinO 
' J 

k j s i 

n0cos0  / 

1 (',->14) 

k.si n0cos0 
J 

• 

1 u 

c '' 

0 

2.1  "■ 

\k .s i n0 

J 

k j^sin0cos0‘^ 

k.si n0cosB 
J 

1 

* 

1 .. 

I U COS0 

+ ® ) 

1 

2,Jt“ 

ikjSinB 

\i  sin0 

i kj  s i nO/ 

'4 

1 • 

o 

1* 

1 * ^ 

/ U COS0 

+ ° ) 

1 ' 

2,£ 

ik.  s i n0 

J 

'l  ^ 

\ i s i n0 

i k.  s i n0/ 
J 

4 

where  the  * indicates  the  complex  conjugate,  and, 


" » kzl 


J dk 


kz,  -iks-  -iks.  -iks^  -''>^5ql 
e Le  ^ -he  ^ -e  ^ 4e 
k-k. 


® kz( 

J dk^ 

o 


(B-1) 

(B-2) 

(B-3) 

(B-4) 


(B-5) 


(b-6) 


B-1 


K-ZO^O 


Integration  of  1^  and  1^ 

It  is  apparent  from  (B-^)  and  (6-6)  that  an  evaluation  of  two 
integrals  is  necessary: 


G 


1 


- J 

o 

_ru 


dk 

k-kj 


-ik(Sj+iz) 

e 


G 


2 


-ik(s j+iz) 
e 


(B-7) 


Transforming  and  G2  hy  translation. 


w “ k - kj 
and  by  rotation 

u ■ iw(S|+iz) 

we  have. 


G,,  -u 


■ i kj  s j +kj  z 
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Thus,  by  applying  the  formula  for  exponential  integral,  we  have 


G,  -e 


G2  “ e 


-ik^Sj+kjZ 


lE,[-ikjS,+k^z]  - i . s,^ 

i2n  1’  *1^ 


-ikj^S,+k^^Z 


(B-8) 

(B-9) 

(B-10) 


(B-11) 


Due  to  the  above  expressions  (B"  1 1) , the  and  1 in  (B-5)  (B-6) 

for  z*0  , are 

4 -ik,s 

1 

m* 


'4  " 


•p  / ,\ii>+l  3n\.  r..  -i  i2n,  >Q 

(-1)  e U,L-.k3S^]  - 0 i • "m  < 


4 


^V|  e **  *"  |E,[-iki^s^J  + j^ni.  ^ 0 


(B-12) 


“m^  l^lt--*Vm^  - ’o'i 


i2TJi 


-ik.s 

i V'  / « \ 1^1  4 m 

- (-1)  e 


8 

m«5 


(B-13) 
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APPENDIX  C 


CALCULATION  OF  INTEGRALS  FOR  COMPLEX  k 


We  evaluate  the  integrals  complex  argu- 

ments In  Eqs.(46  , (51)  and  (52). 

Consider 


G 


3 


I dk 


-ik(s^+iz) 

e 


(C-1) 


where  kj»  complex,  k^*  complex  conjugate  of  k^  and  lmkj<  0.  And  Rekj  ^ 0 


Translating  and  rotating  by 

w ® k-kj  and  u = iw(Sj+i2) 

we  have 

i®(s,+iz) 

I -u 

G,  = r — du 

-ikjSj+kj2  ^ 

For  kj,  with  2“0,  we  have  the  following  relations: 

If  s > 0. 

Re(-ikjS,)  < 0 
Im(-ikjSj)  < 0 

If  s,<  0. 

Re(-ikjSj)  > 0 

Im(-ik^s^)  ^ 0 


for  Re 


kj^  0 


for  Re 


k3>0 


(C-2) 


(C-3) 


Cl 
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For  with  z“0,  we  have  the  following  relations: 

If  s,>  0, 

Re(-ikj^s,)  > 0 

Itn(-ikj^Sj)  > 0 for  Rek^^  ^ 0 
If  s,  < 0, 

Re(-ik^^Sj)  < 0 

Im(-ik^^Sj)  > 0 for  Rek^^  > 0 

In  addition  to  the  above  for  both  k^  and  k^^  with  2*0, 

Re(  i s joo)  * 0 

lni(is,co)  <0  for  Sj  ^ 0 

Thus,  applying  formulas  of  exponential  i ntegral s ^ we  have  the  following 
results  for  z«0. 


For  Re  kj  ^ 0 


s,  ^ 0 


e 3 •|E,[-ik3S,J  - J } . s,  ^ 0 


(C-i4) 


and  for 

Re  k|^  ^ 0 


e ’U,[-ik^^s,]  + 1 . s ^ 0 


-ik.s, 

e U,[-iki^s,]  + Q 


(C-5) 


Thus  for  2=0,  the  1^  and  are  given: 


Re  kj  > 0 


>5- 


44  V 1 s 

0 (-1)"^’  e 3 ^ r_!L  - T i2n 

m=l 


- T i 


s > 0 
m < 


(Cont'd) 
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APPENDIX  D 

EVALUATION  OF  EXPONENTIAL  INTEGRAL  OF  COMPLEX  ARGUMENT 

To  calculate  the  exponential  integral  of  a complex  argument  E,[-ik.s  ] 

I X m 

as  given  in  Eqs.(51)  and  (52),  we  put 


2 “ -i k ,s 
X m 


(D-1) 


where  k^  represents  a complex  wave  number  k^  or  k^^  and  the  s^  has 
been  given  in  Eq.(A-5),  Appendix  A. 

.17 


The  exponential  integral  is  given  by 


EjCz]  « -Y-  log  2 -E  (iarg 


Z - Lj 
n“l 


2I  < n) 


Let  the  above  2 be  given  in  polar  coordinates 
ia 

2 = r e 

, lm(-ik.s  ) 

. . I I • I I . “1  i m 

with  r = -ik.s  , a = tan  •--/  r 

' X m'  ’ Re(-ik^s^) 

Substituting  2 Eq.(D-3)  into  Eq.(D-2),  we  have 

EiC-ik^s^]  = -C(r,Qf)  - iS(r,a) 


where 


/ • * n n 
C(r.»)  . V * I09  r +E  t-'i 

n=l 


® / I \ n n 

S(r.„)  - « + S X-li—I. 

' ' n . r\  • 


no’ 


n»1 


n*n: 


where  V = 0.5772. 

Introducing  a new  variable  3 , 

P * o-  + tt 


(D-2) 


(0-3) 

(D-4) 


(D-5) 


(D-6) 


(D-7) 
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we  have 


z'  = (-1)"  [r"cos  nP  + i r^sin  ng] 
, Hence,  Eq.(D“5)  is  given  by 

] = -C(r,3)  -iS(r,B) 


(0-9) 


where 


C(r.3)  = Y + log  r + 5 Llcosn_B 

n • n • 


S(r.3)  = (6-n)  + £ n B 

n=l 


(D-10) 


The  exponential  integral  of  a complex  argument  -ik^s  ) 

in  Eq.(O-l)  approaches  the  exponential  integral  of  a pure  imaginary  argument 

-'•‘Vm  0 (k^j  ^ 0); 


^i 


fork  >o.s„>0 


(0-11) 


where  Ci  and  Si  represent  the  cosine  and  sine  integral,  respectively. 

Evaluation  of  the  exponential  integral  of  pure  imaginary  argument 
is  carried  out  by  making  use  of  the  standard  subroutine  "SICI"^  which  has 
been  given  in  the  Scientific  Subroutine  Package.  The  program  computes  the 
sine  and  cosine  integrals 

S i (x)  and  C i (x)  , x > 0 

and  thus  by  making  use  of  the  above^  we  compute  the  exponential  integral 
E|(ix)  » -Ci(x)  + i[Si(x)  - , x > 0 

and  Ej(-ix)  with  x > 0 by  the  complex  conjugate  of  Ej(ix). 


D2 
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Evaluation  of  Eqs.(51)  and  (52) 


S i nee 


-ik.s  s 


'1  lit  •^r  m 

e = e e 


i ni 


(D-12) 


we  have 


-ik.s  k£.s 

X m r • 1 T •*'1  m 

e E,L-iksJ=-e 

1 £ m 


j[cos(k  s )C(r,g)4-sin(k,  s )S(r,g)] 
x-f  m X-  rn 


+ i[-sin(k^^s^)C(r,e)+co5(k^^sJS(r,3)  ]1 

(D-13) 


-i2ne 


and 


-ik,s  ko.s 

3 m ^ ^ I ■ nir  . 


= -2Tie  [sin(k-3  s ) + i cos(ko  s )]  for  k-j  >0,  s > 0 

m m Jr  ni 


(D-14) 


-ik.s  k^^.s 

i 2ne  = 2Tie  ' ’’’[sinfki.  s ^+i  cosfki.  s )]  for  ki,  > 0,  s < 0 

m'  4^  ’ m 


= 2Tie  ' '’’[sin(kr  s )+i  cos(ki^  ■ 
4r  m 4^ 


(D-15) 


Substituting  Eqs.(D-l3)  to  (D-15)  in  (51)  and  (52),  one  obtains  the 
expressions  for  l5(kp  - I^Ck^  . 


-ik,s 


5 "iE,[-ik3,J  - 'f  I , J 0 , Re  k,  > 0 


= -e 


k-^.s 

J I m 


|cos(kj^sJC(r,3)  + sin(k2^,sj[s(r,g)+ 


-ie  "’|-sin(k^^sJC(r,3)+  cos  (k^^s  J [S  ( r ,3) | (D-16) 


0 

® ’2tt},  s^  ^ 0 . Re  k^  > 0 


* -e  "’jcos(k^^sJC(r,3)  + s i n(k^^  ^s  J [s  ( r ,3)  - 

-ie  "'|-sin(k4^s^)C(r,3)+  cos(kj^  ^s  J [S  ( r ,3)  - 2n]j 

D3 


(D-17) 


r 


I 
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Since  Ecj.(D-l3)  is  fonnolly  i(l*'ntical  to  liie  complex  conjuyate  of 

(4’or  Appendix  of  Reference  19,  and  since  the  formula  gives 

expressions  for  both  small  and  large  r for  a rapid  numerical  calculation, 

the  expression  of  the  ui^  +i'lo  has  been  used. 

’Op  '^r 


For  small  r,  Eq.(D-13)  remains  unchanged. 

For  large  r.  Eg. (0-13)  takes  the  following  form 


- 1 k ,s 

X m . r ■ I T 
e E,L-'k.s  J 

I X m 


N 

(n- 1)  ! 
n»I 


cos  n b 

n 

r 


N 

+ iD  (n-l)! 
n-l 


s i n n p 


(0-18) 


where  N ■ truncation  number. 

Integration  of  Eqs . 46,  59  and  60 

Integration  with  respect  to  0 in  (46)  (59)  and  (60)  has  been  carried 
out  numerically  by  employing  the  trapezoidal  rule.  For  the  singular  behavior 
which  appears  in  the  neighborhood  of  0 ■ ^ , a change  of  variable,  adopted 
in  Reference  2,  has  been  tsnployed  for  convenience. 

«• 

cos‘^0 


which  gives 

0'  -*  0 for  0 0 

0*  -*  00  for  0 — 

2 

indicating  the  singular  behavior  has  been  shifted  to  6'  “ at  which  the 

Integrand  tends  to  zero. 

For  the  Integration  with  respect  to  O'  , the  upper  limit  has  been 
truncated  at  77  as  being  sufficient. 

In  the  neighborhood  of  the  critical  point  0^  , however,  when  the  t 
values  nearly  satisfy  the  relation  (29),  singular  behavior  occurs  also 
wfwi  k^-kj^-*0as  It  Is  soon  from  Eq.  (26b): 


D4 
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s ■ ^ 


k sec  0 •/ 1 - 4tcos0 
o 


()  - 4t  cosO)  -•  0 for  0 -•  0 


In  these  cases  the  Lagrange  interpolation  method  has  been  employed. 
The  integration  has  been  assumed  to  have  the  form 


I . Jim. 

e -26 

c 


whore  P ■ 0.00033  has  been  finally  adopted  and  M(0) 
of  Eq.  (46)  or  (59)  or  (60)]. 


(0-0  )^.  [ Integrand 
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APPENDIX  E 

FLOW  DIAGRAM  AND  PROGRAM  SYSTEM 
With  the  following  input, 

IN  PROGRAM  IN  TEXT 

TU  T “ reduced  frequency 

i^^f 

fan  |S^  " discharge  rate  for  a pair  of  fans 


SK 

k 

S 

nozzle  flow  factor 

ROW 

Pw 

ei 

water  density 

SELA 

A 

o 

» 

seal  leakage  at  the  equilibrium 

BH 

h 

a 

cushion  height 

DS 

A 

a 

displacement 

i 

XNF 

n 

a 

number  of  pair  of  fans 

1 

SL 

L 

a 

nominal  ship  length 

SB 

B 

- 

nominal  ship  beam 

A 

a 

a 

half  cushion  length 

B 

b 

a 

half  cushion  width 

FH 

l^n 

a 

Froude  number 

one  obtains 

the  following  results: 

RV 

V 

r 

a 

real  part  of  the  volume  coefficient 

i 

IV 

V. 

1 

a 

imaginary  part  of  the  volume  coefficient 

RA 

a 

r 

real  part  of  the  escape  area  coefficient 

lA 

a. 

t 

imaginary  part  of  the  escape  area  coefficien 

i 

AMP 

S/rf 

a 

amplitude  of  the  bubble  pressure  coefficient 

PHA 

phase 

lag 

- 

phase  lag  of  the  bubble  pressure 

{ 

ADM 

^a 

a 

added  mass  coefficient 

DAM 

a 

damping  coefficient 

Appendix  E includes  three  sections;  Flow  Diagram  and  Program  System  for 
Subroutine  VLLl,  INPUT-OUTPUT  SAMPLE  and  PROGRAM  LISTING. 
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APPENDIX  E 
(Cont 'd) 

INPUT  OUTPUT  SAMPLE 


R-2040 


It 


RUN  ACUP 


DEFINITIONS 
RV  - REAL  PART  UOLUME 
IV  = IMAGE  PART  VOLUME 
RA  = REAL  PART  ESCAPE  AREA 
lA  = IMAGE  PART  ESCAPE  AREA 
AMP=  DUBLE  PRESSURE  AMPLITUDE 
PHA-  DUPLE  PRESSURE  PHASE 
ADM=  ADDED  MASS  COEFFICIENI 
DAM=  DAMPING  COEFICIENT 


INPUT  FAN»SKfROW 

»SELAi-BH»DS» 

XNF 

0.03 

0.65  1.94  0.333 

0.72  298.  2. 

INPUT  SLrSB.AfD. 

FN 

13.83 

2.12  6.915  1.06 

0.72 

>> 

INPUT  TUf  IF  ANY 

. OTHERWISE 

KEY  IN 

-1 . 

10. 

TU  = 

10.00  RV  = 

-0.073 

TV 

0 . 00  1 

RA  = 

-0.064 

lA  ■•= 

0 . 003 

AMP= 

0.439 

PHA^ 

-150.130 

ADM  = 

-.2013E-02 

DAM  = 

2177E-01 

» INPUT  TUf  IF  ANY.  OTHERWISE  KEY  IN  -i. 

-1 . 


END  OF  EXECUTION 


V. 
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00100 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

01000 

01100 

01200 

01300 

01400 

01500 

01600 

01700 

01800 

01900 


03100 

03200 

03300 

/ 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04800 

04900 

05000 

05100 

05200 

05300 

05400 


CCCC  MAIN  PROGRAM  FOR  COMPUTION  OF  OOLUM  FLUX  AND 

CCCC  ESCAPE  AREA  OF  ACO  HEAUING  IN  UNIFORM  STREAM 

CCCC  IN  ADDITION  CALCULATE  BUDLE  PRESSURE^  ADDED  MASS 
CCCC  AND  DAMPING  COEFFICIENTS 

cccccccccccccccccccecccccccccccccccccccccccccccccccccccccccccccc 

PROGRAM  COMP 

COMMON  /01/SL»SBr AfB»FN7TU»EPSI >GAM»BTA 
COMMON  /02/BETf ZTl »ZT2rZT3»ZT4 

COMMON  703/  ALFl  r ALF2  f ALF3 » ALF4  » NST12  H'ISTIS  » ILn 
COMMON  /04/NTSl fNTS2rNTS3f KTl»KT2,KT3  7KTi2»K  fi3 
COMMON  /Cl/  EPS(3) 7GAB(3) 

DATA  EPS/. 001  f .0001 7 . 0001/ 7 GAB/ . 03325 7 .011.7 .011/ 

DATA  BTA7ZTI7ZT27ZT37ZT4/.O57 .047.067.172./ 

DATA  ALFl 7 ALF2  7 NTS 1 7 NTS2  7 NTS3  7 KTl / . 033  7 . 05  7 20  7 20  7 30  7 10/ 
DATA  GAMA7PA7ROA7G/I.472048. 70.00237732.2/ 

ALF3=ZT3 

ALF4=ZT4 

KT2=NTS2 

KT3=NTS3 

CCCC  SL  = NOMINAL  SHIP  LENGTH 

CCCC  SB  = NOMINAL  SHIP  BEAM 

CCCC  A = HALF  CUSHION  LENGTH 

CCCC  B = HALF  CUSION  WIDTH 

CCCC  FN  = FROUD  NUMBER 

CCCC  FAN  = FAN  DISCHARGE  RATE  FOR  PAIR  OF  FANS 

CCCC  XNF  = NUMBER  OF  PAIR  OF  FANS 

CCCC  SELA  = SEAL  LEAKAGE  AT  THE  EQUILIBRIUM 

CCCC  BH  = BUBLE  HEIGHT 

CCCC  ROW  = WATER  DENSITY 

CCCC  D3  = DISPLACEMENT 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
TYPE  600 

600  FORMAK'  » DEFINITIONS'/'  RU  = REAL  PART  VOLUME' 


02000 

CCCC 

SL  = 

02100 

CCCC 

SB  = 

02200 

CCCC 

A = 

02300 

CCCC 

B = 

02400 

CCCC 

FN  = 

02500 

CCCC 

FAN 

02600 

CCCC 

XNF 

02700 

CCCC 

SELA 

02800 

CCCC 

BH  = 

02900 

CCCC 

ROW 

03000 

CCCC 

DS  = 

FORMAK 


IV  = IMAGE  PART  VOLUME'/ 

RA  = REAL  PART  ESCAPE  AREA'/ 
lA  = IMAGE  PART  ESCAPE  AREA'/ 
AMP=  BUBLE  PRESSURE  AMPLITUDE'/ 
PHA=  BUBLE  PRESSURE  PHASE'/ 


ADM=  ADDED  MAi 


COEFFICIENT'/ 


TYPE  300 
FORMAT ( ' 
ACCEPT  100 7 
TYPE  400 
FORMAK  ' 
ACCEPT  IOO7 
BL=2.*A 
BB=2.*B 
CONTINUE 
TYPE  500 
FORMAK  '• 
ACCEPT  100 7 
IFdU.LT.O. 
FORMAT (7F) 


DAM=  DAMPING  COEFICIENT'/) 

» INPUT  FAN73K7ROW7SELA7 BH7DS7XNF' ) 
FAN  7 SK  7 ROW  7 SELA  7 BH  7 DS  7 XNF 

» INPUT  SL73B7A7B7FN' ) 

SL7SB7A7B7FN 


» INPUT  TU7  if  any.  otherwise  KEY  IN  -1.') 
TU 

) CALL  EXIT 
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05500 

05600 

05700 

05800 

05900 

06000 

06100 

06200 

06300 

06400 

06500 

06600 

06700 

06800 

06900 

07000 

07100 

07200 

07300 

07400 

07500 

07600 

07700 

07800 

07900 

08000 

08100 

08200 

03300 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

10000 


OR^O . 

01=0, 

CALL  OLLKOR^OIf  I ) 

01  = 0. 

AR=0 . 

CALL  OLLI  (AI  f AR>2) 

IFdU.EQ.O.)  GOTO  20 
SP^FN^SORKBL^G) 

SlG-TUiJ^G/SP 

SIG2--=S1G#SIG 

FC=SK/SIG 

BL2=BL.'!<BL 

CA=BL*BB 

OOL^CA'KBH 

PO-DS/CA 

P00(3=PQ/0QL 

FCTOR=PO/BH-'i(CA 

COF-= . 5*R0W*BL2*BL 

C0P=.5*RaW*SP*&L2 

Cl=GAIiA#<  1 .iPA/PO) 

C2=SaRT(2.)kP0/R0A) 

C3=SaRT(2.;J:pO>KROA) 

C4=BL2/R0W/G 

C5=BL/R0W/G 

aR=0R*C4+FC*AI*C5*C2 

ai=0I;)tC4-FC>KAR*C5*C2tFC*SELA/C3+XNF>XFAN/3IG 

DLQR=POOO;KaR 

BL(:M=POOO*ai 

PB=  ( 1 , -C 1 *DLaR  ) *’«2+  < C 1 *DLQI ) ^^2 
PBR=-Cl!K<  1 .-C1*BLQR)/P& 

PB I =-C  1 ^C 1 -KBl-Q  I /PB 
ADMC=PBR*FCT0R/SIG2/C0F 
DAMC=--PBI*FCTOR/SIG/COP 
AMP=SQRT  ( PBR*PBR  fPBI>f:PB.I ) 

PHAS=ATAN2  ( PB  I r PBR ) >!<S7 , 29577 
TYPE  200  » TU » OR , 01  f AR  f AI  f AhP  f PhlAS  » ADMC  f DAilC 
200  F0RMAT(2Xf 'TU  = 'F8.2»'  RO  = 'F10.3f'  10  =‘'F10.3/ 

1 15X»'  RA  = 'F10.3f'  lA  ='F10,3/ 

2 15Xr'  AiMP=  'F10.3»'  PHA=  -FIO.G/ 

3 i5Xf'  AtiM=  'E10.4,'  DAH=  'E10.4/) 

GOTO  10 

20  AhP=0. 

PHAS=-90. 

TYPE  200 » TU » OR  f 0 1 r AR  f A I » AMP » PHAS 

GOTO  10 

ENB 


E9 
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00100 

CCCC 

SUBPROGRAri  FOR  COhPUTATION  OF 

00200 

CCCC 

VOLUME  FLUX»  IF  IP  = 1 

00300 

CCCC 

ESCAPE  APEAf  IF  IP  = 2 

00400 

SUBROUTINE  VLLI < ZC » ZS f IP ) 

00500 

DIMENSION  THR(3) F IHL ( 3 ) f FRR ( 3 ) f 1 

00600 

COMMON/MYTEX/  SM(H> 

00700 

DIMENSION  ZCC(2) fZ33(2) fFTR(3) fI 

00800 

COMMON/LOG/Cl fC2fR1 fR2 

00900 

COMMON  /V1/3Lf3BfAfBfFNf  IUfEPSI 

01000 

COMMON  /V2/BETFZTiFZT2FZT3FZT4 

01100 

COMMON  /V3/ALF1 f ALF2  f ALF3  f ALF4  f i 

01200 

COMMON  /V4/NTS1 fNTS2fNTS5fKT1fK 

01300 

COMMON  /CFl/  TVfSI fS2fS3fS4fSGC: 

01400 

COMMON  /CFll/  3K05fTTMfTM 

01500 

COMMON  /CF2/  EP ( 4 ) f El ( 4 ) f EPD ( 4 ) 

01600 

COMMON  /Cl/  EPGS<3) fGAB(3) 

01700 

DATA  X11fX44fX15/11. f144. f315./ 

01800 

ZCC(1)=0. 

01900 

ZCC(2)=0. 

02000 

ZSS(2)=0. 

02100 

ZSS(1)=0. 

02200 

P=3. 14159 

02300 

TP-=2 . *P 

02400 

P2=P*P 

02500 

SLZ^SLHtSL 

02600 

SB3L=  • 5^3E</SL 

02700 

ABL2=4.YA*B/SL2 

02800 

G=32.2 

02900 

3LMN=77. 

03000 

BLF=0. 000333 

03100 

EPSI=EPSS( IP) 

03200 

GAM=GAB( IP) 

03300 

BET=GAM 

03400 

LP=3 

03500 

LP1=LP+1 

03600 

X=-A 

03700 

Y=B-EPSI 

03300 

GAM3=3.*GAM 

03900 

KT12=KT1+KT2 

04000 

KT13=KT12+KT3 

04100 

RMA=SLMN-BET-ZTl.'lt(KTl-l)-ZT2*KT 

04200 

M00=RMA/ZT4 

04300 

RMB=^M00*ZT4 

04400 

DELTA=RMA-RMB 

04500 

KT14=KT13+M00 

04600 

KT15=KT14H 

04700 

NT312=NTSlfNTS2 

04800 

NTS13=NTS12+NT33 

04900 

IF(TU.LT. .25)  GO  TO  127 

05000 

GC=1 ./<4.*TU) 

05100 

TCR=ACOS(GC) 

05200 

TVCR  =SIN ( TCR ) /COS ( TCR ) **2 

05300 

EPS=LP*BLF 

05400 

RTA=SLMN-TVCR-EPS-ALF1*(NTS1-1 ) 

05500 

L0=RTA/ALF4 

-MLr2.tNT:32-ALF;5KH  rS3 


ElO 
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05600 

05700 

05800 

05900 

06000 

06100 

06200 

06300 

06400 

06500 

06600 

06700 

06800 

06900 

07000 

07100 

07200 

07300 

07400 

07500 

07600  127 

07700 
07800 
07900 
08000 
08100 
08200 
08300 
08400 
08500 
08600 
08700 
08800 
08900  C 
T 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

10000  2216 
10100 
10200 
10300 
10400 
10500 
10600 
10700 
10800 
10900 
11000 
11100 


RrH::^L0i.ALF4 

DELTH-=RfA--RT& 

NTS14^NTS13fL0 

Nr315=NTS14fl 

RA=^T0CR-EPS-GAM3 

LQH^RA/BTA 

R£i=LQH*BTA 

DEL==RA-RD 

LaMi=LaMfl 

GHA-^TOCR-BET-ZTl^CKTl-i  )-2T2i^KT2-ZT3t-KT3 

LOO=SHA/ZT4 

3ME<-LOO*Zr4 

DELTC  =3MA-3nB- , 005 

LT14=KTi  31-LOO 

LT15=LT14fl 

IT3T^2 

Kril5“KTi5 

IF(TVCR.GT.77»  ) ITST-=1 
IF(T0CR,LE.66.  ) I TST==0 
IF(IT3T.EQ,2)  KTli5=LT15 
CONTINUE 

SFC:=1 ./( 180.*BLF.t*2) 

SFCC=1  ./<180.*0AM*.K2) 

SP=FN*SQRr(G*SL) 

3IG™TUKG/SP 

SIG2=3IG)«*2 

3P2-^3P>K.K2 

FN2=FN**2 

SGC2=G/SP2 

FCT^=4 . >!<P2.'I(3P21:SL2 

IF (IP. EG. 2)  FCT=FCT/SL 

FCTR=FNK>i(4/(2.*P2) 

IF  ( IP.  EG.  2)  FCTR=2.*P2M<3L 

'IT'  STANDS  FOR  CONTR  REPET  OF  INTGR  FOR  X--X: GIVES  RSUL 

DO  1 IT--=1»2 
CT3--0. 

CT4=0. 

SII=0. 

SIR=0. 

ZCl^O. 

ZS1'=0. 

DO  2 NST=lf6 

IF(NST.EO. 1 .0R.NST.EQ.6)  GOTO  2216 

IF(ITST.GE.l)  GOTO  2 

CONTINUE 

IF(NST.EG.l)  TV=BET 
IF(NST.EG.2)  TV=&LF-fTVCR 
IF ( NST . EG . 3 ) TV  =TVCRFEPS+BET 
IF(NST.EQ.4)  TV=TVCR-EPS 
IF (NST. EG. 5)  TV=GAM 
I F ( NST . EG . 6 ) TV=TVTDEL-GAri 
IF(NST.EG.6.AND. ITST.GE. 1 ) TV=GAH 
IF (NST. EG. 2)  THR(1)=TV 
IF (NST. EG. 4)  THL(1)=TV 
IF(NST.EG.1>  ILM=KT15 


Ell 


f 

t 

R-2040 

11200 

IF(NST.EQ.2.0R.NSr.EQ. 4)  ILH=3 

11300 

IF(NST.EQ.3)  ILM-NTS15 

11400 

II- <N3T.Ea.5)  ILM=^3 

11500 

IF  ( N3  f . EO . 6 . AND . I CST . EO . 0 ) ILii^LOiil 

11600 

IF(NST.EQ.6. AND. ITST.GE. 1 ) ILH^KTii5 

11700 

DO  3 I=lfILH 

: 

11800 

IN==I  + 1 

. 

11900 

CALL  COEFKTU) 

12000 

IF(IP.EQ.2)  CSN-C3N/SGC2/CN 

12100 

IF(IP.E0.2>  SGCO=SGCO.KCN 

12200 

FDrSG-=FDr*SGCO 

12300  122  CONTINUE 

• 12400 

321==S2-S1 

12500 

543=34-33 

12600 

R1=X;KCN  + Y*3N 

12700 

R2=X*CN-Y*3N 

{ ‘ 

12800 

Cl=A*CN  + B.tc3N 

12900 

C2=A*CN-E()1:3N 

13000 

3M( 1 )=R1-C1 

13100 

3M(2>=R1-C2 

13200 

SM(3)=R1+C1 

r 

13300 

SH(4)=R1 fC2 

13400 

SM(5)=R2-C1 

1 

13500 

Srt(6)=R2-C2 

13600 

3M(7)=R2+C1 

13700 

3M(3)=R2+C2 

13800 

CALL  ALO ( I P . 0 i » 02  T 03  f 04  f 05  r 06 ) 

13900 

IF(TU.Ea.O.)  GO  TO  100 

14000 

51S=31.*.'!<2 

14100 

523  = 52*!K2 

14200 

333  = S3;!c.'{<2 

14300 

34S  = 341cH<2 

14400 

310=313*31 

14500 

320=323*32 

1 4600 

330=333*33 

1 

14700 

340=343*34 

i 

14800 

F3=3IG2/( 31*32) 

i 

14900 

FC33=F3 

15000 

F4= ( SP*S 1 *CN-SIG ) **2/ ( -3 10*321 ) 

15100 

F5= ( SP*32*CN-3 IG ) **2/ ( -320*321 > 

; 

15200 

Fi=-3GC0*(SIG~3P*31*CN) **2/ (313*321 ) 

t ; 

[ I 

15300 

F2=3GC0*  < 3 IG-3P*32*CN ) **2/ ( 323*32 1 ) 

15400 

IF( rU.LE. .25)  GOTO  1231 

' t 

15500 

IF ( N3T . GE . 4 . OR . N3T . EQ . 1 . OR . 343 .LI.. lE-7 ) GOTO  123 

15600  1231  FC3'=SIG2/(33*S4) 

[ i 

15700 

FC4= ( SP*33*CN+SIG ) **2/ ( -330*3  43 ) 

! 

15800 

FC5=  < 5P*S4*CN+3IG ) **2/ ( -340*343 ) 

15900 

IF ( IP . EQ . 2 ) F3=SGC0* ( 3IG  + 3P*33*CN ) **2/ ( 333*343 ) 

16000 

IF ( I P . EQ . 2 ) F4=-SGC0* (31 G f 3P*34*CN ) **2/ ( 3 43*343 ) 

16100 

IF( IP.EQ.2)  FC4=-3IG2*3GC0/( 33*34) 

f 

16200  123  CONTINUE 

16300 

IF ( IP . EQ . 2 ) FC3=S IG2*SGC0/ (31*32 ) 

' 

16400 

IF(NST.Ea.2.0R.N3T.EQ.3)  GOTO  19 

16500 

IF(N3T.GE.4)  GOTO  33 

1 

\ 

16600 

CALL  YTEXKSlfl.EPR^EPI) 

16700 

CALL  YTEX  1 ( 32 . 2 r EPR » EP I ) 

16800 

IFdU.LE.  .25)  GOTO  66 

i 

• ! 

EI2  - 

R-2040 


.16900 

CilILAC-lTI.rA 

17000 

CALL  C(jr.F2(l»2) 

17100 

tF(IP.Ea.2)  (3010  t33 

17200 

CALL  VFa(F3fF4tF5»Q1  »Q2»iri»0'!  »u5rU6. 1 »2rYCl » Y31 ) 

17300 

GOTO  305 

17400 

133 

CONTIriUE 

17500 

YC 1 ^F  H<  ( EP  < 1 ) +EPD 1 > > iF2 1 < EP  ( 2 > t-EPD < 2 > ) 

17600 

Y31=--F1  ii(EI<  1 ) lE.iU';  1 > ) fF2i<(EI  (2)  {•E1D<2>  ) 

17700 

Y C 1 = Y C 1 - ( F 1 Y F 2 > -X  Q 1 f F C 5 r U 3 

17800 

YS 1 ==YS1  - ( F 1 PF2 > .Xa2+FC5.Xa  » 

17900 

YCl--=YCl*Fur 

18000 

YS1--=Y81.XFDT 

laioo 

GOTO  305 

18200 

19 

CALL  YrEX2<S3r3rEPRfEPI ) 

13300 

CALL  YTEX2  ( 84  ..  4 r EPFi , EPI  > 

18400 

CALL  C0EF2(3f4) 

18500 

IF(rp.EQ.2)  (3010  191 

18600 

CALL  VFX  ( FC3  F FC 4 F FC5  f Q 1 r U2  f Q3  f u 4 f 05  f Q6  f 3 f •(  f YC  1 f Y31  ) 

18700 

GOTO  192 

18800 

191 

CONTINUE 

18900 

YC  1 ^F3.{;  ( EP  ( 3 ) TEPD  ( 3 ) ) +F4-X  ( EP  ( 4 ) f EPO  < 4 > ) 

19000 

YS1==F3I((EI  (3)  }E1D(3) ) fF4-X(EI<4>f£IIi(4) ) 

19100 

YC7  =-  (F3+F4  ) KOI  } FC4>X03 

19200 

YS7'=  ( F3f  F4  ) ;XQ2-FC  4.YG4 

19300 

YCl=(YCl  + YC7);XFDf 

19400 

YSl  -=(  Y31TYS7)  Y.FEiT 

19500 

192 

CONTINUE 

19600 

IFCNST.EO.S)  GOTO  506 

19700 

FRR  ( I ) ■■=  YC  1 .t  < THR  < I ) - TUCR ) 

19800 

F I R ( I ) = YS  1 < THR  ( I > -TOCR ) m3 

19900 

ALF=BLF 

20000 

GOTO  306 

20100 

506 

CONTINUE 

20200 

DELAC=DELTB 

20300 

GOTO  305 

20400 

33 

TM~- TH 

20500 

CTU=TM.X*0.5 

20600 

S3R=8K05;XTTM 

20700 

S4I~SK05*CTU 

20800 

34R=33R 

20900 

S3I~-34J 

21000 

IF(IP*EQ.2)  GOTO  335 

21100 

CALL  KMOL ( S4R  f S4 1 f 3 1 G f SP  f CN  f FD 1 f FD2  f FD3  f FD4R  f F&4 1 > 

21200  • 

GOTO  336 

2130C 

335 

CALL  KMP(S4RfS4I fSGCQfSIGfSPfCNfFKI fFK2fACfBD) 

21400 

336 

CONTINUE 

21500 

CALL  YTEX3(S3RfS3If3fEPRfEPI> 

21600 

CALL  YTEX3  < 34R  f 34 1 f 4 f EPR  f EPI > 

21700 

CALL  C0EF2<3f4) 

21800 

IF(IP.EQ.2)  GOTO  337 

21900 

YC3--=FD  4R*  ( EP  ( 3 ) +EPD  ( 3 ) TEP  ( 4 > f EPD  ( 4 ) ) -FD4 1 ♦ ( E 1 < 3 ) 

22000 

1 fEID<3)-EI(4)-EID(4)) 

22100 

YS3-Fri4R*  ( El  ( 3 ) +E ID(  3 ) +EI  < 4 ) +EID  ( 4 ) ) +FD4I * ( EP  < 3 ) 

22200 

1 +EPD(3>-EP<4)-EPD(4) ) 

22300 

YC3  = < FD  1 *Q  1 f FD2.Xa3  f FD3.tQ5  ) 

22400 

Y58~  ( -FD  I *02-FD2;X04-FD3-XQ6  ; 

22500 

YCi-  ( YC3  1 YCO  ) -XFDTSG 

r 


R-2040 


JL’AOO 

f 5 1 V iSOfYOO  ) *Fli  1 .G 

22/00 

GO  III  220 

22000 

22/ 

YC2  ^Al'f  (LI'v  2 ) ILfli-.  2 ) ) -Ulit  v t 1 v • 1 1 1 IH2  . ) ( 

22000 

1 or,  * U 1 M.  4 ) I L'l  lu  4 ) > 1 1 1 ( <1  ) 1 1. 1 11 1 .1  / > 

22000 

Y52  O<liM[0‘l2)  HOlu*)  ) lOG  KU,I  i W HMluO)  > 

22100 

1 - on  ^ u.  !•  ( 1 ) 1 FPu  14))  foi.  1:  ( r 1 V 4 ) h:  uu  4 . ) 

22200 

•jpo  -1  K 1 .101  II  KJ'KO  5 

22200 

Y50  1 l\  1 YOO  -f  K2Y04 

22‘>00 

Yi'.l  Ol‘:2l'ri;0).Yl  HI 

221.00 

Y'.l  V i02l)fS0)H  HI 

22*00 

220 

roN  1 1 oiji: 

22/00 

II  <oor.Gr.5)  OHIO  500 

22000 

1 Kl  \ 1 ) Yi;  l ,v  ( 1 III . ( l ) - 1 Oriv ) i |;2 

22/00 

1 11  1 l ) - YGl  M 1IIL(  l )--l  VCIO  rlcO 

2 ->000 

01  1 t:|  1 

2U00 

00  10  ^v)* 

/•»200 

500 

CON  1 INOi: 

2-K500 

II  V [ ISl  .C-L.  1 ) 00  VO  125 

2-»  400 

11  uN;')  1 . 1.0  . .') ) OllVO  12* 

24000 

1 OKI  1 ) 0 I NI  rO)Y(t2 

2 4*00 

I 1 IvU  )~YC  U 1 lllv'i  1 1 

24  700 

1 I 1 i 1 ) ^ YSlKlMKli  1 ) 

•24000 

Oil  OOH 

2 4'/ 00 

00 ro  20* 

2’jOOO 

126 

CONI  INOI  : 

25100 

HIILOC  HIO 

25200 

OHIO  205 

25200 

125 

CON  1 INOr 

25400 

II  \ MSI  .1  0.  1 ) HLI.OO-  HI  1...  1 0 

25500 

11^  1 lor.l  11.2)  HCLOC-HEL  1 C 

2 5/.  00 

00  10  20:. 

25/00 

i: 

COOIMIIOI  ION  OF  O-FI.X  FOO  Ills.  25  IN  HLLlHJ 

25000 

66 

CON  1 INOI:. 

2 5'?  00 

COI  L Y VI:.\2iC>2.2.FFK'.FFI  ^ 

2*000 

COI  L Y ii:x2v54, 4.1  i-k.fim:  ) 

2A100 

coi  l.  COI  021  1.4) 

26200 

HFl.Oir^HLL.  10 

2*200 

iri[F'.L0.2)  GO  10  *65 

26  400 

YC  1 -^F-IXFI'l  1 ) FEPHi  1 ) )-F5'Na:pi  2)  H;ptu2)  ) 

2*500 

YO  1 F4ia:  1 V 1 ) IFlln  1 ) )--F5fiL  1 i 2) il  i:rn2)  ) 

2**00 

YC2-iFC  4 i ( EP  ( 2 ) FEPH  i 2 ) ) -FC5*  ( EP  < 4 < ♦ l:.PH  v 4 ; ) 

26/00 

Y02=  f F C 4 ♦ i E 1 1 2 ) ♦ E 1 H <3'>>  -FC5  i E n 4 ) I E 1 H i 4 ' 

2*000 

YC  1 -YC  1 ♦ ( F4  IF5)'V0l  1 1-F4  .VO  1 1F5 1 52  ) 10)2102^05 

2*900 

Y5  1 =Y5  1 1 i -14  FF  5 ) '102  1 i -F4'15  1 FF5I(52  > 104  I F 2 10* 

27000 

Y C 2 Y C 2 1-  ( - F C 4 F F C 5 ) '1 0 1 F ( - F C 4 1 5 2 1 1 ■ C 5 r 5 4 ) 1 0 . 5 1 F C , o 5 

27100 

Y S 2 - Y S 2 - i - F C 4 f F C 5 ; 1 0 2 - < - F C 4 1 5 2 1 F C 5 1 5 4 > 1 0 4 - 12:  * + 0 * 

27200 

YCl  -iY(:iFYC2)*l  DVBG 

27200 

Y5l^(  YSl  FYB2)'lFltTSG 

27400 

GO  10  205 

27500 

6*5 

CONI t NOE 

27600 

YC  1 -1  1 >1  ( EP  ( 1 ) FEPH  ( 1 ) ) iF21  ( EP  ( 2 ) l EPH  ( 2 ) ) 

27700 

Y;3  1 =F  1 >1  i E I i 1 ) 1 E [ D ( 1 ) > IF  2 1 1 E t ^ 2 ) 1 E 1 D l 2 ) ) 

27000 

YC2'rFJia  P(2)  lEPDv  2Y  ) FF 4 1 1 EP  ( 4 ) FEPD  ( 4 ) ) 

2/900 

Y52  02* <E 1 (2) FEIIK  2) ) IF  41 lE 1 ( 4) FEllK  4 ) ) 

20000 

YC  1 :»YC  1 - 1 F 1 FI  2 ) 10  1 FFC2  102 

20100 

Y5l  ’T!U“iFl  FF2)102FFC2.104 

20200 

YC2'>'YC2  ■(F2FF  4)101  FFC4102 

EI4 
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28300 

28400 

28500 

28600 

28700  100 

28000  C 

23900 

27000 

29100 

29200 

29300 

29400 

29500  305 

29600 

29700 

29000 

29900 

30000 

30100 

30200 

30300 

30400  306 

30500 
30600 
30700  3 

30800  991 

30900 


Y52-=  rS.l  1 ( F3  IF  4 > .)'L'.'-FC4  t.U4 
YCl-^(  IC  I fYC2)  .l«Fii  T 
YSl“(YSlfY'32)  ltH  Ul 
GOrO  305 
CQMT  INUi: 

COL  OF  ULL  FUR  T0U--0. 

CtU.L  Y IFX 1 ( S2 » 2 ? I :RR-  » FP I ) 

COLL  CULF2(..S2) 

YC  l-CGN'KFC  l t<  M'C  l M lU  ( FP  v 2 ) lin  UU  ; -GJPU.')  -U  I ) 

I F ( I f • . C 0 . 2 ) Y G I - C C)  N F 0 V .'K  ( E I ( 2 ) ♦ F ;n  U . > ) 1 R.’ ; 
ir([p.E0.2)  YCl^O. 

ItFl  OC-=liELrO 
CUN  riWUE 

[ F ( I . EU  . I . UR  . I . |-:u  . 1 LM  ) VC  1 ^0 . 5 f(  YC  1 
rF< I .EU. 1 ,UR. I.FU.ILH)  YGl-0.5iYSl 
[F(l.CU.l)  OLF-KEl 

II  ( [ .Ell.  1 .ONG.  I rSf  .EO.O.ONU.NG)  .EU.6>  OLF--iiLL 
ZC I ZCl  lYC IFOLF 
ZGl  -ZGl fYGIKOLF 
Nsn-Nor 

OLF-OLFC i IN » NS  fl > DELOC  y 1 I ST ) 

TV--=iO*OLF 

[F  i NS T . ECl . 2 ) ri  IR  I I f 1 ) = I 0 
IFiNSf  .Ell.  4)  TIILi  1 11  ) =10 
CUNflNLIE 

FUKh0Tl2Xy2E10. 4) 

IFi I U.LL.O. 25)  Guru  11 


31000  1 1 k NS  I .EU  . 5 )C  ('3  ■Sl'CC'^l-Xl  1 tF  I R'.  ) 1X4  4 4^1"  1 K v 2 >-A  151F  i R i I ; ; 

31100  IF  (NS  r .Ell.5)CT4"SFr,C4<<-Xl  I'M  1 [ ( 3)  1X4  4.'M-  T I i2;-X  15Y1  il  i 1 > ; 

31200  ir(NSf.NE.6)  GOTO  2 

31300  S1R1=X  imvFRLl  1 )-FRR(3)  ) 

31400  3iR2--X444ilFRL(2)-FR'R(2)  ) 


31500 

31600 

31700 

31800 

31900 

32000 

32100  2 

32200  11 

32300 

32400 

32500 

32600 


S 1 R3  = X 1 5 Y ( FRL ( 3 ) -FRR ( 1 ' ) 

S I R = S 1 C Y ( S I R t f S I R'  2 1 3 1 R 3 ) 

SI  [1--=X1  ItlFlLi  1 )-FlR(3>  ) 

SI  l2=-X44'lc<FIL(2)-FlR(2)  ) 

St  I3^^-X  l5t(FIL(3)  -FIRi  1 ) ) 

S 1 1 =SFC4(  ( S 1 1 1 +3 1 1 2 IS  1 1 3 ) 

CONTINUE 

ZC  = ( ZC I f SIR 1CT3 ) /FCT 
Z3~(2S1 ISII  1CT4)/FCT 

IFtTU.EQ.O.  .AND.  [p.Ell.2)  ZS^ZSl/TClR 
X=-X 

ZCC(  IT)=--ZC 


32700 

ZSSdD^ZS 

32800 

IF  (IP.  Ell.  2)  GOrO  99 

32900 

1 

CONTINUE 

33000 

99 

CONTINUE 

33100 

ZC  - ( zee  ( 2 ) -zee  ( 1 ) ) ;!(2  . - 

33200 

ZG=(Z33(2) -ZSS( 1 > ) *2. 

33300 

1 F ( I P . E 0 . 2 ) Z S :=  Z S - S D 1 i L 

33400 

IFi  TP.  El).  2)  Ze^=~Ze-ADL 

33500 

RE  TURN 

33600 

END 

E15 
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00100 

SUt<l'.‘UU  TIME  KHF'  ( ? S 4 1 r SbCG r 3 1 G !> SP » CrJ » FK 1 r FK j 7 

00200 

1 ACftiD) 

00300 

3 

F0l<iiAr(6E12.4) 

00400 

XX^2..J?3P.'K3Ii3.KCN 

00500 

3If32^=S:[G.K*2 

00600 

S4l<2=--S4R.'i<  K2 

00700 

S4;[2  = !3  4I.}cjK2 

00800 

XY2==S4R2F34I2 

00900 

XP=2,.^<S4R 

01000 

XY4^XY2*;!c2 

01100 

F K 2 ^ - 3 (3  C □ « S I G 2 / X Y 2 

01200 

AA  ( SP*S4R;^<C^^^S :t G ) .''..K2-  ( SP:{<34 1 .YCi'J ) **2 

01300 

BB-2  , t ( SP-'i<34R;};CN  f 3 1' G ) Y3P+.3 4 1 

01400 

DD-^2.;{<34I.'«(S4R2-S412) 

01500 

CC=^4  7.t34R:K34I2 

01600 

C:CDD=CC.i<*2}DD*.!<2 

01700 

A C - ( A A Y G C - B B ■!<  B D ) / C C B B * 3 (3  G G 

01800 

BB=="-(BBtCG+AA*BB)/GGBB*3GC0 

01900 

FK1:^=2..'I<AC 

02000 

RETURN 

02100 

ENB 

02200 

SUBROU  r I NE  ALG ( I P 7 Q 1 7 02  7 03  7 0 4 7 05  7 06  > 

02300 

GOrtMON/MYTEX/  Sri<e) 

02  400 

BIMENS.ION  ASM(3>  7AL(8)  7GI  (3) 

02500 

COMHON/LOG/C 1 7 C2  7 R 1 7 R2 

02600 

C 

01-=RE(  .11  HIB)  702-'=;[n<  IlfllB)  703^-=RE(  :i;3f;[3B>  704=^Inri3Fi3B) 

02700 

c 

05^^4^;E  ( 191-1 9B  ) 7 06^  1 h ( 1 9 T 1 9B ) 

02300 

p;[=3. 14159 

02900 

BG  1 M^-=^lr8 

03000 

ASH  ( ii ) ~AB3  ( SH  ( M ) ) 

03100 

IF(ASH(h))  37374 

03200 

4 

AL(H)^ALGG(ASM(M) ) 

03300 

GG  ro  1 

03400 

3 

GDMTIrlUE 

03500 

1 

GGNTINUE 

03600 

aP  1 = -AL  ( 1 ) + AL  ( 2 ) -AL  ( 3 ) -f  AL  ( 4 ) 

03700 

AP2---AL  ( 5 ) -AL  < 6 ) fAL  < 7 ) -AL  ( 3 ) 

03800 

AP=AP1 fAP2 

03900 

AP3=AL ( 3 ) -AL ( 1 ) 

04000 

AP4=AL(2)-AL( 4) 

04100 

AP5^AL<5)-AL<7) 

04200 

AP6-AL<3)-AL<6) 

04300 

BG  31  K=^l7  8 

04400 

IF<SH(K) .LT .0. ) Gi(K)=PI 

04500 

IF(SM(K) .GE.O. ) G1(K)=0. 

04600 

31 

CONTINUE 

04700 

AOl-^^Gl  (3)-Gl  ( 1 ) 

04800 

Aa2=Gl(2)-Gl (4) 

04900 

AQ3  = G1  <5)-(31  (7) 

05000 

A04^G1 (8)-Gl (6) 

05100 

3II1-~~G1  ( 1 )+Gl  (2)-Gl<3)+Gl  ( 4> 

05200 

SI  I1B=:G1  (5)-Gl  (6)  fGl  (7)  -31(8) 

05300 

02'=ST  t I f 31 1 ID 

05400 

S 1 1 3=R  1.  * AP  1 -C  l;»!AP3  ~G2*AP4 

05500 

SI  I3D=R2.tAP2-GH<AP5-C2*APo 

E16 
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05600 
05700 
05800 
05900 
06000 
06100 
06200 
06J00 
06'»00 
06500 
06600 
06700 
06800 
) ) 

06900 

07000 

07100 

07200 

07300 

07400 

07500 

07600 

07700 

07800 

07900 

08000 

08100 

08200 

08300 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 


f}4  SI  [31811311 
Ul-iM' 

EU<l3  = -Fvl  1G[  11  [C  Mail  I 1C2mHI2 
8l<I.^lt=-^<2KS  1 1 ILHCI  K0a3iC.’*(Ul4 
U 3 ~ S R 1 3 + R 1 3 l.i 

if([p.i;q.2)  K't  ri.iRii 

SR9-- . 5 XR'l  .'i'*2  .'tt'iPl  'R1 1C  1 ♦ wM,  ( 1 j -i  ■'  3 > ) t R i li,.M  v .U.  v 2 ; -ill.  ( 0 ^ 

) 10 . 5 *C  1 1 12*  V nL  (1)1  nl.  \ 3 > > . '.j  i l-..‘  •-  (lL  ^ 2 ) t nL  v -t ) ; 

S19  = - . 5*Rl.l;.ii2.*5i  1 1 -Ri  tCl  li  U1  O ) -ul  ^ 6 ; 

IR  I *C21n  G I < 2 ) -(31  ( 4 ) ) 10 . 5 IL  1 1 *2  1 M)  1 i ) ( G I ^ 3 ; > - . ;/tC2'1  f 
(G1  (2)  IG  I ^ •'( ) ) 

r.R9V.i  = -() . 5 1R2>  12  li  11-2  llv2*i:  I * ^ A1  v 5 > -nl.  k / > ) {<2 1<i;2  ^ \ GL 
(6)-AL(0)  ) -0.5'|(Cl*!!2t^GH5)  lOL  <,  7 ; ) 10 . 5 *i;2*  12 1 ^ i iH  6 ,•  I nl.  ^ 8 

S 1 9 li  ■ - . 5 * R 2 * * 2 1 Sill  i 1 1 R 2 K L I lA  G H 5 ) --  n U 7 > ) - R 2 1 C 2 1 u.i  U 6 ) 

-G I ( 8 ) > -0 . 5*C 1 * 12*  uJ  1 ( 5 ) lu l ( 7 ) > 1. SIC 2 **2*  u3 1 ^ 6 / I u 1 Ui  > ; 

G5=:3R91SR9D 

06=^3 1 9 IS  1 9Li 

RETURN 

lINti 

SUtiROU  r L NE  RrtOL  ( S4R  » 3 4 I » 3 IG  * SP » CN  mVl  r A2  r ('iS  r f) 4 !>  f'llj  > 
COhPUlEG  CGEF  OF  KOHPLEX  WAVE  NUAHERS 
A 1 ■ = F m > A 2 F D 2 » r)  3 = F D 3 » A <>  I • D 4 R » A 5 F U 4 .[ 

31G2-S [G**2 
SP2-SPMK2 
CN2-^CN.i<*2 
X~G4R 
Y=:34  [ 

X2-^A'112 

X3^X2*X 

Y2-Y**2 

Y3'=Y2*Y 

Y4-Y3KY 

XY2-X21Y2 

XY2P-XY2i<*2 

XY3P--XY2113 

XZ=XY2*2. *X 

XY=3.  .1X2-Y2 

A-^3P*3IG.tCN 

H=3P2*CN2 

U--6.*X2*Y212.*Y4 

0=2,*Y*X3-6.*X*Y3 

U~Ulc*2 10**2 

G - ( SP*X*CN 13G ) **2 - ( SP*Y  *CN ) **2 

H =2 . * < SP*X*CN  IS ,IG ) *3P* Y*CN 

A3^-31G2/XY2 

A4=(G*U+H*0)/W 

A5-(G*0-H*U)/U 

A2=-2.*(A4*XiA5*Y) 

A1^-2.*A4 

RETURN 

END 

SUDROUT I NE  YTEX 1 ( 3K  » K » EPR . EP 1 ) 

COMMON/M YTEX/  3H ( 8 ) 

D I MENS I ON  EPR ( 4 » 8 ) » EP 1 1 4 . 8 ) 
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1 iOOO 

12 

1 t JKi  iA  r \ 2X  » ' C . 3 F 1- 1 F S I '■  / ) 

1 1100 

13 

1 oi-.n.ii  woi:i2.-)) 

1 1200 

1 4 

rOKhnI  U^IXf  'SKiI'/  > 

11300 

r^3. 14 159 

ll'lOO 

n -4  K2  . 

1 1^.00 

in)  n-lFil 

1 1600 

3I\H  iMifi  ( SiK  !<Sn  i rt  ) > 

1 1 -■’00 

C CDS  <3l\n) 

lUtOO 

S SlN\l)KM> 

1 1900 

coil.  S 1 1. 1 ( S 1 » C 1 F SKrt  F 

12000 

If-  ( ;>n  \ 11 ) ) F 3 • 9 

12100 

7 

EI-K  (KfM)=:-C'XC1“S*S1  - I I ts 

12200 

r:p[U\FM)--s*ci-ci'oi  - 

12300 

gd  i d * 

12‘)00 

3 

CDNI INDL 

12500 

GD  ID  6 

12600 

9 

Lf'K'dx  Fii)-^-  -c  m:  l -StS  I 

12200 

tPl  IK  Fh)--S.l'Cl  ICtSl 

12B00 

6 

CONI INUC 

12900 

11 

PDIshO)  (/3Xf  'SINGDLOK'/) 

13000 

RC 1 URN 

13100 

END 

13200 

SDHRD1.H  INi:  YriT\2lSKFRFi:PRFt.Pl  ) 

13300 

C D ri  m J I'l  1 \ Y T i:  X / S 11  \ 0 ) 

13'U)0 

n [ hPNS  I ON  f.'PR  1 4 F 3 ; F I P U ) F 3 ) 

13500 

12 

lilK'i-l.')  1 ( / 2X  F ' C F S F C 1 F SI  ' / ) 

ISaOO 

13 

PDRDOVi  IOC  12. 4> 

13700 

1 4 

FiJRNn  1 K / 2X  F ' SKn  ' ) 

13300 

P 3. 14159 

13900 

IP  ^2.  M' 

14000 

i n . ) 6 il  ■ 1 F 3 

14100 

SKD  -i'iOS  V SK^Sh  >>  n ) ) 

14200 

C-CDS(SI\ii) 

14300 

S "S 1 N V SKil ) 

14400 

COLL  SICI iSI fCI fSKN) 

1 4500 

lF<Sh(M))  7f3  f 9 

14*00 

7 

El ) 'R-  U\  F h ) =:  - C K C 1 - S % S 1 

14700 

t F'  1 ^ K F i'i ) ~ - S K C I 1 C S 1 

14300 

IP (K. CO, 4)  GO  TO  19 

1 4900 

GO  TO  6 

15000 

3 

CON  1 1 nue: 

15100 

GO  1 0 6 

15200 

19 

LF'R  < K • M > ^CPR  \ K F i1 ) -S*  TP 

15300 

Cl  - 1 1 R F M ) -t P r ( R F n > 1 C i I P 

15400 

GO  ro  6 

15500 

9 

EPRiRFM>=-C')tr:i-S'FSI 

15600 

EPl  (RfH)- St-CI-CtSI 

15700 

[F(R.E0.3)  GO  TO  13 

15800 

GO  ro  6 

15900 

18 

LPR  < R F H ) =^F:PR  < R f M ) - I PtS 

16000 

LP I ( R F rl  > -EP  I<  R F H ) -iP*C 

16100 

6 

CONI  [nue: 

16200 

1 1 

porno r</2XF 'SINGL'/) 

16300 

RL  rURN 

16400 

ENG 
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i 


V. 


latiOO 
I A£>00 
16700 
16^00  2 
16900 
I 7000 
17100 
17200 
17300 
I ,M00 
I 7500 
17600 
17700 
171J00 
I 7900 
IJIOOO 

laioo 
10200 
10300 
10  100 
10500 
10600 
10700 
10000 
10900 
19000 
19100 
19200 
I ‘>300 
19-100 

19500  20 

19600  2b 
19  ’00 
19000 
19900 
20000 
20100 
20200 
20300 

20400  30 

20500 
20600 
20700 

20800  31 

20900  26 

21000 
21100 
21200 
21300 
21400 

21500  70 

21600 
21700 
21000 
21900 


1 1 1 Mb'  Y I bXJ  \ 'IIJ » 0 I » K » r.hK  , tPl  ; 

Li  IMLMUIOr)  Ll  Ki  4 fO  ? rbPI  ( -1 .0) 

Cl'.iihlli'l/h4  rtX/  OrtvO) 
rUlvhAl  ( iX»  l0UI0.‘l  ) 

F'-3.  I 4 1 59 

ri-rn2. 

II  0 

It  ( K . I il . 4 . ONli . SR  . IJ  r . 0 . ) I r-3 
IF  ( K . im  . 3 . AND  . SR  . 0 I . 0 . > I I'-- 1 

no  60  h^ifO 

IF  V Sh  V A ) . Fit . 0 . ) 00  to  60 
X — ‘SRtSM^M) 

XY^-X 

Y-S  I FSiKM) 

X2^X»t2 

Y2^'t*F2 

AX-AHSvX) 

E<A-'^X2FY2)  **0.5 
C^OOV.vAX) 

S=3tNv AX' 

lF(Y.OE.OO.)  Y=HO. 

IF (Y. Lb. -00. ) Y^-00. 

E-t;XI'<  Y) 

HLi  = «A 

Ar  -ATilW2uX.  Y) 

I t F 

IFrSK.LQ.O.  ' nC-f 
HL.-0. 

nn^o. 

l F V,  L<A-6  . ) 20 » 30 » 30 
rtri  i 

LU  l<Li  F nil  *l.OS  1. 1 LOA  I <,  iMn  ) YL<C  > 
liL  ICI  niiFS  IN  (FLOA  I (.NhltnC) 
liLi-li[!*nA*FLOA  r ^ hh ) /FLOA  f i,  unn  t I ) in2  > 

MM  MM  F I 

IFdifi.Or.O.  1E-D*HA)  OO  to  2b 
nn - ( 0 . 577  I ALOG  ( IIA ) fnii ) 
tlt:-^llEFA^) 

GO  m 26 

no  .H  liM  - 1 . 5 

nil  ■ nil  1 003  { FLOA  I t hh  ) *80  ) /nil 
nO-iiE  FS  IN  I FLOA  F ^ Ah ) *nO  ) /HH 
nn^nA-Fun/FLOArirtfi) 

COM  I INUE 
CONt INUE 

IFxSR.EO.O. ) GO  ro  69 

GAS*  rF*E*S 

GAG* TF* Etc 

IF(liA-6.)  70.77.77 

IF(FiA-6.)  ro.  77.  77 

GA -L  * < 0 * ( - nn ) f S I GN ( S . X Y > * ( -HE ) ) 

On=E*(-S[GN(‘S.XY)tv-Hn>  IG*v-nE)  ) 

If  V K . Lll . 3 . AND  . SR . G I . 0 . « ANH . Sri  i n t .01  . 0 . ) i iA ' OA—i.n tS 
I F ( K . FO . .1 . ANH  . SR  . G 1 . 0 . . ANH . mu  ri ; . u ) . 0 . ; 0;<  Gii  ~GAi> 
I F V K . FO  . 4 . ANH . SR  . G 1 . 0 . . ilNH . SM  v n ; . L I . O . ) i iii  "( lA-Gi  i:.> 


Fn 


) 


J 
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22000 

I F ( K . EQ  . 4 . AND  . SR . i > I . 0 . . AND . Sri  ( ri  > . L 1 . 0 ♦ > 

22100 

GO  TO  00 

22200 

77 

IF( 1 r .EO, 1 .ANu.GM(H) .GT.O. ) GO  TO  92 

22300 

IF< lT.EQ.3.AND.Sh(rt) .LT,0. ) GO  TO  93 

22400 

GO  TO  94 

22500 

92 

GA^-BD-E.nP*S 

22600 

GB-DE-E+.TP^C 

22700 

GO  TO  80 

22800 

93 

GA  = -Dli-E;{<TP>KS 

22900 

GB=^DE+ETTP'i<C 

23000 

GO  TO  ao 

23100 

94 

CA=^-DD 

23200 

GB-DE 

23300 

GO  TO  30 

23400 

69 

CONTINUE 

23500 

GA=^-E)TBD 

23600 

GB-=0. 

23700 

80 

EPR  ( K , M ) =^GA 

23800 

EPt<KH’'i)=GB 

23900 

60 

CONTINUE 

24000 

99 

FORMAT ( /2X » ' 5INFULAR ' / ) 

24100 

RETURN 

24200 

END 

1 


